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Abstract 

We describe a new universality class for unitary invariant random matrix ensem- 
bles. It arises in the double scaling limit of ensembles of random n x n Hermitian 
matrices Z~\\ det M\ 2a e~ NTrV( - M UM with a > -1/2, where the factor | detM| 2a 
induces critical eigenvalue behavior near the origin. Under the assumption that the 
limiting mean eigenvalue density associated with V is regular, and that the origin 
is a right endpoint of its support, we compute the limiting eigenvalue correlation 
kernel in the double scaling limit as n, N — > oo such that n 2 ^(n/N — 1) = 0(1). 
We use the Deift-Zhou steepest descent method for the Riemann-Hilbert problem 
for polynomials on the line orthogonal with respect to the weight \x\ 2a e~ NV ( x \ 
Our main attention is on the construction of a local parametrix near the origin by 
means of the ^-functions associated with a distinguished solution of the Painleve 
XXXIV equation. This solution is related to a particular solution of the Painleve 
II equation, which however is different from the usual Hastings-McLeod solution. 
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1 Introduction and statement of results 
1.1 Unitary random matrix models 

For n G N, N > 0, and a > —1/2, we consider the unitary random matrix ensemble 

Z^ N \detM\ 2a e- NrFlV ^ dM, (1.1) 
on the space A4(n) of n x n Hermitian matrices M, where V is real analytic and satisfies 

lim - ,\', = +oo. 1.2 

x^±oo \og( X 2 + 1) V ' 

This is a unitary random matrix ensemble in the sense that it is invariant under conju- 
gation, M i — ► UMU^ 1 , by unitary matrices {7. As is well-known [Til [38] , it induces the 
following probability density on the n eigenvalues x\, . . . , x n of M 

n 

P^ N \x u ...,x n ) = Z^ N Y[\ Xj \ 2a e- NV ^ J] \ Xl - x,! 2 . (1.3) 

j=l i<j 

The eigenvalue distribution is determinantal with kernel K n> ^ built out of the polynomials 
Pj,N(x) = Kj^x^ + Kj t N > 0, orthonormal with respect to the weight \x\ 2a e~ NV ^ xS> 
on M.. Indeed, as shown by Dyson, Gaudin, and Mehta, see e.g. [H] EH 38J, for any 
m = 1, . . . ,n — 1, the m-point correlation function 

I POO POO 

R^ N \ Xl ,...,x m )= / ■■■ / P^ N \x 1 ,...,x n )dx m+1 ---dx n (1.4) 

[n m). > /_ 00 J -co 

is given by 

i^'^i, ...,x m ) = det (i^fe ^•)) 1 < iii < m , (1.5) 

where 

n-1 
3=0 

In the limit tt,, TV — > oo, n/TV — ► 1, the global eigenvalue regime is determined by V as 
follows. The equilibrium measure \xy for V is the unique minimizer of 

WilA — I I l°g i rdfi(x)dfi(y) + / V(a:)<i//(x) (1.7) 

JJ \x-y\ J 

taken over all Borel probability measures [i on M. Since V is real analytic we have that 
\iy is supported on a finite union of disjoint intervals [13], and it has a density /v such 
that 

lim -K nN (x,x) = p v (x), x^O. 

n,N^oo,n/N^l 71 

The limiting mean eigenvalue density is independent of a. 
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The factor | detM\ 2a changes the local eigenvalue behavior near 0. This is reflected 
in the local scaling limits of K n ^ around that do depend on a. If is in the bulk of 
the spectrum and pv(0) > 0, then instead of the usual sine kernel we get a Bessel kernel 
depending on a [37]. If is in the bulk and py(0) = p' v (0) = 0, p' v (0) > 0, then the local 
scaling limits of the kernel near are associated with the Hastings-McLeod solution of 
the Painleve II equation q" = sq + 2q 3 — a [S] . 

In this paper we study the effect of a in case is an endpoint of the spectrum which 
is such that the density py vanishes like a square root at 0. For a = the scaling limit is 
the well-known Airy kernel, see the papers [21 ESI ESI 112] and also [31 [T2], and so we are 
asking the question: What is the a-generalization of the Airy kernel? 

For a > —1/2, we have found a new one-parameter family of limiting kernels as stated 
in Theorem 11.11 below. In Theorem 11.11 we also assume that the eigenvalue density py is 
regular, which means the following. 

• The function i h 2 Jlog|a; — s\py(s)ds — V(x) defined for x G K, assumes its 
maximum value only on the support of py. 

• The density py is positive on the interior of its support. 

• The density py vanishes like a square root at each of the endpoints of its support. 



Theorem 1.1 For every a > —1/2, there exists a one-parameter family of kernels K^ d9e (x, y; s) 
such that the following holds. Let V be a real analytic external field on R such that its 
mean limiting eigenvalue density py is regular. Suppose that is a right endpoint of the 
support of p v so that for some constant c\ = c\y > 

Pv{x) ~ — M 1//2 as x — ► — . (1.8) 

7T 

Then there exists a second constant ci = C2y > such that 

lim T \ r ,K n , N ( T A-j^, r J ^) = K^{ X) y; s) (1.9) 
whenever n, N — > oo such that 

lim n 2/3 ( — - l) = L e R (1.10) 
and s = — c 2 yL. 



For a = 0, the limiting kernels reduce to the kernel 

K ed9e (x y s) = Ai (z + s)Ai'(y + s) - Ai'(x + s)Ai(y + s) 



x-y 



which is the (shifted) Airy kernel from random matrix theory mentioned above, see also 
Subsection 14.11 below. For a^fl, a new type of special functions is needed to describe 
the limiting kernel K^ d9e (x, y; s). This description is given in the next subsections. 
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Figure 1: Contour for the model RH problem. 



1.2 The model RH problem 

We describe K^ d9e (x,y; s) through the solution of a special Riemann-Hilbert (RH) prob- 
lem, that we will refer to as the model RH problem. 

The model RH problem is posed on a contour E in an auxiliary £-plane, consisting of 
four rays E x = {arg£ = 0}, E 2 = {arg£ = 27r/3}, E 3 = {arg£ = tc}, and E 4 = {arg£ = 
— 27r/3} with orientation as shown in Figured! As usual in RH problems, the orientation 
defines a + and a — side on each part of the contour, where the +- side is on the left when 
traversing the contour according to its orientation. For a function / on C \ E, we use f± 
to denote its limiting values on E taken from the ±-side, provided such limiting values 
exist. The contour E divides the complex plane into four sectors Qj also shown in the 
figure. 

The model RH problem reads as follows. 



Riemann-Hilbert problem for \I/ a 

(a) ty a : C \ E -> C 2x2 is analytic. 

(b) * Q>+ (0 = *a,-(C)fj |Y force E x , 
* Q)+ (C) = * a ,-(C) (Li ^,forCeS 2 , 

¥^+(0 = *a,-(C) J) , for C G E 3 , 

(c) * Q (C) = C" CT3/4 75 f ] (I + 0(l/C 1/2 ))e-(i^ /2+ < 1/2 ) CT3 as C -> oo. Here a 3 = 
( o -i ) i s t ne third Pauli matrix. 

(d) * a (C) = O ( \ C X a \ C X a ) as C - 0, if a < 0; and 
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*auH ) L | ,A ifa>0. 



icr icr 

Here, and in what follows, the O-terms are taken entrywise. Note that the RH problem 
depends on a parameter s through the asymptotic condition at infinity. If we want to 
emphasize the dependence on s we will write ty a ((] s) instead of ^ a (C)- 

The model RH problem is not uniquely solvable. Indeed, if *$? a is a solution, then 
()7 1) is also a solution for any 77 = T)(s), and it turns out that this is the only freedom 
we have (see Proposition 12.11) . 

Theorem 1.2 The model RH problem is solvable for every sGl Let be a solution 
of the model RH problem and put 



(1.12) 

y Q ,+ (x; s)e- ama * [ Jl, /or x < 0. 
T/ien t/ie limiting kernel K^ d9e (x,y; s) can be written in the "integrable form" 

Kedge{ a) = Mx-,s)My-,s)-M*-,s)My-,s)_ 

2ni{x — y) 

The function ip 2 depends on the particular choice of solution \1/ Q to the model RH 
problem. Indeed, for any 77 we have that the mapping ty a 1— > ( * ? ) leaves ^1 invariant 
and changes ^2 to ^2 + However, this does not change the expression (I1.13P for the 
kernel K%f 9e (x, y; s). 

It follows from (I1.12p and part (c) of the model RH problem that ip\ and ip2 have the 
asymptotic behavior 

iM*; s) = - 7 =L__ e -S- 3/2 -^ 1/2 (l + 0(x^/ 2 )), (1.14) 
^(*; s) = !^ e -f^ 2 — 1/2 (i + O^- 1 / 2 )), (1.15) 
as x — ► +00, and 

s) = v^lxr^cos (^|x| 3/2 - s|x| 1/2 - avr - tt/4^) + 0(x~ 3/4 ), (1.16) 

ip 2 (x;s) 
as x — ► —00. 



V^xl^sin Qlrcl^-slxl^-oTr-Tr/^ + 0(aT 1 / 4 ), (1.17) 



Remark 1.3 The kernel K^ d9e (x,y; s) describes an edge effect for the random matrix 
ensemble ( 11.11) . If we assume that is the rightmost point in the support of py, and if given 
M we let A max (M) denote its largest eigenvalue, then it follows under the assumptions of 
Theorem II. 1[ in particular the limit assumption (11.101) . that 

lim P ((cm) 2 / 3 A max <t)= det (l - K a>s \ (too) ) , (1.18) 
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where Ka,s\( too \ is the trace class operator in L 2 (t, oo) with kernel K ed9e [x, y; s). To 
prove (11.181) one must show that the operator with kernel ^ c ^ 2/3 K n ^ ^ ^^2/3 ? ( Cl „)2/3 j 

converges in the trace class norm on L 2 (t, 00) to the operator with kernel K edge [x, y; s). 
This requires good estimates on the rate of convergence in (11.91) . which can be established 
as in [T2] . 

For a = 0, the kernel is the (shifted) Airy kernel, and the Fredholm determinant (11.181) 
has an equivalent expression in terms of a special solution of the Painleve II equation. 
The resulting distribution is the famous Tracy- Widom distribution [321 H3] . It would be 
very interesting to find an analogous expression for general a. The connection to the 
model RH problem given in Theorem 11.21 can be used in obtaining such an expression, 
following the approach of [5] and [27] . We are planning to address this question in a future 
publication. 

1.3 Connection with the Painleve XXXIV equation 

The model RH problem is related to a special solution of the equation number XXXIV 
from the list of Painleve and Gambier [29], 

„<W + 2 s „ + M 2 -< 2a ) 2 . (1,9, 

2u 

All solutions of (j 1.19ft are meromorphic in the complex plane. 
Theorem 1.4 Let ty a ((;s) be a solution of the model RH problem. Then 

c(^(C; S )e(^""^" 2 )"^(i !: 7)c /4 



u(s) = i— lim 

2 as c^oo 



12J 



1.20) 



exists and satisfies (11.19ft . The function (11.201) is a global solution of ( 11.191) (i.e., it does 
not have poles on the real line), and it does not depend on the particular solution of the 
model RH problem. 

The connection with the Painleve XXXIV equation leads to the following characteri- 
zation of ipx and ip2- 

Theorem 1.5 Let u be the solution of Painleve XXXIV given by (11.201) . Then there 
exists a solution of the model RH problem so that the functions ipi and ^2 defined by 
( 11.12ft satisfy the following system of linear differential equations 



d fip 1 (x;s)\ _ ( u'/{2x) i — iu/x\ fipi(x;s) 

dx~ \Mx; s)) = \-i{x + s + u + {{u'f - (2a) 2 )/(Aux)) -u'/{2x) J \ip 2 (x; s) 



(1.21) 



and have asymptotics ( 11.141) - 01.17ft . 

In fact we will prove that for ( G C\E, 

— * (C- s) - ( M ' /(2C) * ~ m/C \ V (C- s) (1 22) 
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from which (ll.2ip readily follows in view of (11.121) . We emphasize that (jl.2ip and 
(ll.22p hold for one particular solution of the model RH problem. Any other solution 
(77 0) ^a(C; s ) a l so satisfies a system of linear differential equations, but with matrix 

i*7(2C) i-iu/(\(l 0\ 

-<(C + a + tt + ((u') 2 -(2a) 2 )/(4iiC)) -u'/(2C)J\-V I J ' 1 ' 

In order to make Theorem 11.51 a genuine, i.e., independent of the ^ a RH problem, 
characterization of ip\ and ip2, we need an independent of formula fll.20p characterization 
of the solution u(s) of equation fll,19p . This can be achieved by indicating the asymptotic 
behavior of u(s) as s — > 00, cf. the characterization of the Hastings-McLeod solution of 
Painleve II equation [28]. We discuss this issue in detail in the last section of the paper, 
see in particular Proposition 14. II and the end of Remark 14.21 where the possible asymptotic 
characterizations of the solution u(s) are given. 



I. 4 Overview of the rest of the paper 

In Section [2] we give the proofs of Theorem 11.11 and Theorem 11.21 We start by presenting 
the RH problem for orthogonal polynomials on the line [23]. The eigenvalue correlation 
kernel K n ^ can be explicitly expressed in terms of the solution of this RH problem [TTlfTB] . 
As in earlier papers, see e.g. [21 IH El El EEH EE], we apply the Deift-Zhou steepest descent 
method for RH problems, see [IB] . For the local analysis near 0, we need the model RH 
problem for \l/ a (C; s) as introduced in Subsection [L2J We show, following the methodology 
of [25], that the model RH problem has a solution for every s£K. Then we follow the 
usual steps in the steepest descent analysis for RH problems, which lead us to the proofs 
of Theorem 11.11 and 11.21 

Section [3] is devoted to the proofs of Theorem 11.41 and Theorem 11.51 We start by 
discussing the RH problem, in the form due to Flaschka and Newell, associated with the 
Painleve II equation q" = sq + 2q 3 — v. Following [2] , we show that for a special choice of 
monodromy data the Flaschka-Newell RH problem is related to the model RH problem. 
The parameters in the Painleve equations are related by v = 2a + 1/2. The monodromy 
data corresponds to a solution of Painleve II which is different from the Hastings-McLeod 
solution that has appeared more often in random matrix theory [H EJ [9l [28l H2]. The 
known results (asymptotics, Lax pair etc.) for the RH problem for Painleve II are then 
transferred to the model RH problem, and then used to complete the proofs of Theorems 

II. 41 and 11.51 In particular it gives rise to the special solution u of the Painleve XXXIV 
equation defined by fll.20p . 

In Section H] we make some concluding remarks. For the important special cases a = 
and a = 1, we show how the model RH problem can be explicitly solved in terms of 
Airy functions, and how the limiting kernel K^ d9e (x,y; s) as well as the special Painleve 
XXXIV solution u can be explicitly computed in both cases. Our final remarks concern 
the characterization of u, in the case of general a, through its asymptotic behavior at 
infinity. 
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2 Proof of Theorem ED and \T7Z\ 

2.1 The Riemann-Hilbert problem for orthogonal polynomials 

The RH problem for orthogonal polynomials on the line, for our particular weight, is the 
following (cf. [23]). 

Riemann-Hilbert problem for Y 

• Y : C \ R -> C 2x2 is analytic. 

A \x\ 2a e- NV{x) \ 

• Y + (x) = Y-(x) ( g i j for x G E\{0}, with R oriented from left to right. 



Y(z) = (I + 0(l/z))[ Q z _ n ) asz^oc. 



• If a < 0, then Y(z) = O (J j J« J as z -> 0. If a > 0, then Y(z) = O ( \ \ ) as z -»• 0. 
The RH problem has the unique solution 



/ 1 , x 1 /• p ra , j y( g )| g | 2 "e-^W , \ 



Yiz) 



K n,N 27TZK ni jv S — Z 

9 ■ / \ f p n -iAs)\s\ 2a e~ NV{s) , 

> -2niK n _ 1:N p n _ ljN (z) -K n _i,iv / as/ 

\ ./to S - Z / 



(2.1) 



where Pj t j^(x) = Kj^x^ + • • • is the orthonormal polynomial with respect to the weight 
\x\ 2a e~ NV( - x \ By (11.61) and the Christoffel-Darboux formula for orthogonal polynomials, 
we have 

K nN ( x ,y) = \x\ a \y\ a c-* Nmx)+v(v)) Kw ~ 1,JV Pn ' N ^ Pn ~ 1 ^ Z , ( 2 .2) 

«n,iV x - y 

Thus, using (12. ip and the fact that dety = 1, we may express the eigenvalue correlation 
kernel directly in terms of Y: 

K n , N {x,y) = —i-^larl-lyl-e-^W+vft,)) ( x ) Y -\y)Y + {x) Q . (2.3) 

The main idea for the proof of Theorems 11.11 and 11.21 is to apply the powerful steepest 
descent analysis for RH problems of Deift and Zhou [16] to the RH problem satisfied by 
Y . In the case at hand it consists of constructing a sequence of invertible transformations 
Y i— > T i— > S I— > i?, where the matrix-valued function i? is close to the identity. By 
unfolding the above transformations asymptotics for Y and thus, in view of (12. 3p . for 
K n) N in various regimes may be derived. Our main attention will be devoted to the local 
behavior of Y near 0. Around we construct a local parametrix with the help of the 
model RH problem, which we next discuss in more detail. 
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2.2 The model RH problem 

The model RH problem is not uniquely solvable. 

Proposition 2.1 Let ty a be a solution of the model RH problem. Then the following 
hold. 

(a) det ty a = 1. 

(b) For any r\ G K (which may depend on s), we have that (^?) ^ a also solves the 
model RH problem. 

(c) Any two solutions are related as in part (b), i.e., if^fa and are two solutions 
of the model RH problem, then = (J?) ^i 1 ^ for some rj = rj(s) . 

Proof, (a) We have that det \P a is analytic in C \ {0}, since all jump matrices have 
determinant one. In case a < we get from condition (d) of the RH problem that 
det ^(C) — C(|C| 2a ) as ( — > 0. Since 2a > —1 it follows that the singularity at the 
origin is removable. In case a > we find from condition (d) of the RH problem that 
det^ a (C) = 0(1) as ( — > in fij U fi 4 . Thus the singularity at the origin cannot be 
a pole. Since det ty a = 0(\(\~ 2a ) as ( — > 0, it cannot be an essential singularity either 
and therefore the singularity at the origin is removable also in this case. Hence det ty a is 
entire. From condition (c) of the RH problem it follows that det ^ a {C) —> 1 as £ — ► oo, 
and so part (a) of the proposition follows from Liouville's theorem. 

(b) It is clear that ?) \I/ Q satisfies the conditions (a), (b), and (d) of the model RH 
problem. To establish (c) it is enough to observe that 

v? O^TlO =c<T3/4 7lG ( J + 2^ ( i i)) 

= r* /4 ^=() i)a + o(i/c 1/2 )) 

as ( — > oo. 

(c) In view of part (a) we know that ^fa is invertible. Then ^a\^a ) _1 is analytic 
in C \ {0} and, by arguments similar to those in the proof of part (a), the singularity at 
the origin is removable. As ( — > oo we get from condition (c) of the model RH problem 
that 

The statement now follows from Liouville's theorem. □ 

In the following we will need more information about the behavior at the origin of 
functions satisfying properties (a), (b), and (d) of the model RH problem. The following 
result is similar to Proposition 2.3 in [9]. 

Proposition 2.2 Let \1/ satisfy conditions (a), (b) ; and (d) of the RH problem for 
Then, with all branches being principal, the following hold. 
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If a — | ^ No, there exist an analytic matrix-valued function E and constant matrices 
Aj such that 

*(C) = £(C)C aCT3 ^-, forCeQj. (2.4) 
Letting Vj denote the jump matrix for \& on Hj, we have 

A, = A 4 v x , A x = A 2 v 2 , A 3 = A 4 v 4 , (2.5) 

and 

1 1 

A 2 = ( 2cosa7r 2cosct7r I . (2.6) 

• If a — \ G N , then ^ has logarithmic behavior at the origin: There exist an analytic 
matrix-valued function E and constant matrices Aj such that 

9{C) = E{Q * A,, forCettj. (2.7) 

\o r Q / 

Letting Vj denote the jump matrix for \P on E 3 -, we now have 

At = Ai vi, A t = A 2 v 2 , A 3 = A 4 u 4 , (2.8) 

and 

( e 3m ^ 

A 2 = \ I . (2.9) 

\ e 71 "^ 4 e ni/4: 

• In all cases it holds that det Aj = 1 and 

(At) 21 = {AA 2l = 0. (2.10) 

Proof. The statement f 1 2 . 1 j) is an immediate consequence of the explicit formulas for 
the Aj's. 

Consider the case a — | ^ N . Define E by f)2.4p . i.e., let 



E(() = f(C)^ 1 r B », for c e n,-, (2.ii) 

with as in (I2.5p . (12.61) . Then E is analytic in C \ S. We now show that E is indeed 
entire. The relations (I2.5P and the condition (b) of the model RH problem imply that E 
is analytic also on Si U £2 U £4. Moreover, on £3 

ez\Qe + (Q = c a3 A 3 v 3 A, 1 c; a - 

Now, by (12. 51) . (12. 61) . and straightforward computation 

A 3 v 3 A' 1 = A 2 v 2 v^ 1 v 4 v 3 A" 1 = e 2a7Tia3 = (: Qff3 Q a3 . (2.12) 
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Thus, E is analytic also on S3, and therefore in C \ {0}. 

We next show that the singularity at is removable. If a < 0, we see from (12.111) and 
the condition (d) of the model RH problem, that as ( — > 



\C\ a \C\ a \ n fl Anf\C\~ a °\ _n A ICI 



1 2a 



Eic) -°{\c\ a \c\ a )°{i icrj-°vi ici 2 V 

so (since 2a > —1) the isolated singularity at is indeed removable. If a > and ( — ► 
in Qi we find in the same way (also using (^1)21 = 0) that 

bfk i)°(r ,cr)-o; 

so that is bounded near in and thus cannot be a pole. Since cannot be an 
essential singularity either, we conclude that the singularity is indeed removable. 

In case a — ~ G No the proof is almost identical, only now the equation (12.121) is 
replaced by 



-1 -2i 
-1 



A 3 v 3 A 2 1 = A 2 v 2 v 1 1 v 4 v 3 A 

c° y v° c 



C a -h( a logt\ f( a ±C a logC\ [:J)) 



□ 



2.3 Existence of solution to the model RH problem 

We will need that for s G R the model RH problem indeed has a solution. To prove 
existence of a solution to the model RH problem it suffices to prove existence of a (unique) 
solution ty£ pec ' to the special RH problem obtained when the asymptotics (c) at infinity 
is replaced by the following stronger condition 

y(s P ec) {0 = {I + 0(1/C))C-- /4 ^ Q [) e-(K 3/2 +< 1/2 )- (2.14) 
as ( — > 00. 

A key element in the proof of unique solvability of the RH problem for ^a Pec ^ is the 
following vanishing lemma (cf. [25]). 

Proposition 2.3 (vanishing lemma) Let a > — 1/2, sGK, and put 9(C) = 9(C; s) = 
2^3/2 _|_ s( ai/2^ Suppose that F a satisfies the conditions (a), (b), and (d) in the RH problem 
for ty a but, instead of condition (c), has the following behavior at infinity: 

F a (C) = O(l/0C a3/ ^(\ jV 0(C)ff3 , (2.15) 
as C — > 00. Then F a = 0. 
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Proof. The ideas of the proof are similar in spirit to those in Deift et al. 
be defined as follows: 



G«(C) 



F a {C,)e e ^ 



F a (()e 
[ F a (()e 



-1 

1 ) ' 

1 0\ (0 -1 

e 2am e 29 ^ l) [l 

1 0\ 

e -2a7Ti e 2e(0 i I > 



9(0*3 



for C G fii, 

for C g ft 2 , 

for C G Q 3 , 
for (Gfl 4 . 



Let <X 



(2.16) 



Then G a satisfies the following RH problem. 



Riemann-Hilbert problem for G a 

(a) G a : C \ R -> C 2x2 is analytic. 

(b) G a ,+(C) = G a) _(C)«G a (C) for C e R \ {0}, where 

for C > 0, 

(2.17) 

for C < 0. 




V G a (0 \ / -, . e 2a7ri e 2e + (C) 



(c) G a (C) = 0(C- 3/4 ) as C -> oo. 

(d) G a has the following behavior near the origin: If a < 0, 
G«(C)=o(19« asC-0, (2.18) 



icr icr 



and if a > 0, 



Ol J_« LLJ asC^0,ImC>0, 

gm = < (2.i9) 

O ). a y,_ a ) asC-0,ImC<0. 



The jumps in (b) follow from straightforward computations which uses that #+(C) + 
6-(C) = for ( < 0. The behavior (c) of G a at infinity (uniformly in each sector) follows 
directly from (12351) . fl2TT6|) . and the fact that Re 0(C) < for ( G fi 2 U fi 3 . The behavior 
(I2.18P at the origin is immediate from the condition (d) of the RH problem for F a , and 
so is the behavior fl2TT9l if C -> with (eOiUfl 4 . To prove fl2TT9|) if ( ->• with 
C G Sl 2 U ^3, we need Proposition 12.21 Consider first the case a — | ^ N and £ G f2 2 . 
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Then we have, using (12TB . (IXIj) . f[23j) . and fl2~TUj) 

G„(C) = S.(C)e'«>" ( e .J em) J) (J "o 1 ) 



-0(0^3 



e 



where * denotes an unspecified constant. Using the boundedness of E and 9 at the origin, 
we find (12.191) as ( — * in the sector Q 2 - The case C G Sl 3 is treated similarly. Using 
(12. 7p . (12.81) instead of (12 .4p . (I2.5p . the same argument works in case a. — \ £ N . Note 
that in spite of the logarithmic entry in (12.81) . there are no logarithmic entries in (I2.19p . 
Introduce the auxiliary matrix-valued function 

H a (() = G a (Q (G a (C))*, CeC\E. (2.20) 

Then H a is analytic and 

H a (()=0(C 3/2 ), asC^oo. (2.21) 

From the condition (d) in the RH problem for G a it follows that H a has the following 
behavior near the origin: 



O I \ C \2a W\2a J aS C °> in case « < 0, 

HM={ X a ( 2 - 22 ) 
1 I as ( — > 0, in case a > 0. 



Since a > —1/2, we see from (12.211) and (12.221) that each entry of H a + is integrable over 
the real line, and by Cauchy's theorem and (12.211) 

f H a , + (()d( = 0. (2.23) 
Jm 

That is, by (T2T20]) . 

G a , + (() (G a ,40)* d( = 0. (2.24) 
Adding (I2.24p to its Hermitian conjugate and using (I2.17P we obtain 

G a A0 [«Ga(0 + («G«(C))1 (Ga,-(C))*dC 
+oo 

G«AC)(l o)(Ga,-(C))*dC+ [ G a ,.(() ( 2e 'l m °) (G«,-(C))*dC 

(2.25) 
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Here we also used that 0+(C) — ~~ $-(C) G ^ f° r C < 0> which holds because s is real. 
The identity (12.251) implies that the first column of vanishes identically on R. Thus, 
in view of the form of the jump matrix in (12.171) . the second column of G at+ vanishes 
identically on R as well. It follows that the first column of G a vanishes identically in the 
lower half-plane, and the second column vanishes identically in the upper half-plane. 

To prove that the full matrix G a vanishes identically in both half-planes, we shall use 
a Phragmen-Lindelof type theorem due to Carlson [7J UT]. Define for j = 1, 2, 

' \ (GcMC), forImC<0. 

The conditions of the RH problem for G a yield that both gi and g<i have analytic contin- 
uation across (0, oo) and that they are both solutions of the following scalar RH problem. 

Riemann-Hilbert problem for g 

• g : C \ (— oo, 0] — > C is analytic with jump 

9+tt) = 9-(0 e- 2am e 29 -^, for C G (-oo,0). (2.27) 

• 9(0 = 0(C 3/4 ) as C ^ oo. 

• ^(C)=0(|C|- |a| ) asC^O. 

We are going to prove that this RH problem has only the trivial solution. 
Let g be any solution and define g by 

2(C 2 ), ifReC>0, 
g(() = <J g (^ e -2*™ e -2(^+sO^ if ReC < 0, ImC > 0, (2.28) 

#(C 2 )e 2Q ™e- 2( K 3+sC) , if Re C < 0, Im C < 0. 

The jump property (12.271) ensures that g is analytic across the imaginary axis. 
Now define 

H0=(j^y '^(C 4/3 ), forReC>0, (2.29) 

with (as usual) the principal branches of the fractional powers. Then it can be checked 
that h is analytic in Re£ > 0, bounded for ReC > 0, and satisfies 

|^(C)I < Ce- c|c|4 , if C G zM, 

for some positive constants c and C. Hence, by Carlson's theorem, h = in Re£ > 0. 
Therefore g = 0, and so g\ and g2 are both identically zero. It follows that the full matrix 
G a vanishes identically in both half-planes. 

Thus F a = by ( 12. 16ft . and this completes the proof of the proposition. □ 

We now show how (unique) solvability of the RH problem for \fr^ pec ^ can be deduced 
from the above vanishing lemma. 
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Proposition 2.4 The RH problem for W& p has a unique solution for every s Gi. 

Proof. The idea of the proof is this: Given a solution tya P to the above RH prob- 
lem, we show how to construct a solution m a to a certain normalized RH problem (i.e., 
™>a{C) — > / as £ — > oo) characterized by a jump matrix v on a contour E, and con- 
versely. To prove the proposition it therefore suffices to prove (unique) solvability of the 
normalized RH problem. This can be done by utilizing the basic relationship between 
normalized RH problems and singular integral equations. We recall briefly, in our setting, 
some standard facts regarding this relationship. For further details, and proofs, the reader 
is referred to the papers [T7J [HI [191 E]> an d to the appendix of [32J. 
Let C denote the Cauchy operator 

Ch(() = ^- f^-ds, />GL 2 (E), (GC\S, (2.30) 

Z7TZ Jy; s — C, 

and denote by C±h((), £ G E, the limits of Ch((') as — > C from the (±)-side of E. The 
operators C± are bounded on £ 2 (S). Let 

v(C) = (MO)-V(C), CeE, (2.31) 

be a pointwise factorization of t>(£) with v±(Q G GL(2,C), and define u± through 

v±(0=I±u ± ((), CeE. (2.32) 

Our choice of factorization will imply that u± G £ 2 (E) fl L°°(E). The singular integral 
operator : £ 2 (S) — > £ 2 (E), defined by 

CJx = C+(hwJ) + C-(hw+), h G L 2 (E), (2.33) 
is then bounded on L 2 (E). Moreover, it makes sense to study the singular integral equation 

{l-C u )n = I (2.34) 
for fi G / + L 2 (E). For if we write \i — I + h, then (I2.34p takes the form 

(l-C u )h = C u IeL 2 (E). (2.35) 
Suppose that /i G / + £ 2 (S) is a solution of (12.341) . Then, indeed 

m a (C) = 1 + C{fi{u + + w_))(C), C e C \ E, (2.36) 

solves the normalized RH problem. Thus, if we can prove that the operator 1 — C^ is 
a bijection in L 2 (E), then solvability of the RH problem for m a , and hence of that for 
^a Pec \ has been established. Bijectivity of 1 — C u in £ 2 (S) is proved in two steps. We 
first show that, for an appropriate choice of u — (ou-,u+) in the above factorization, 
1 — is Fredholm in I/ 2 (E) with index 0. Second, we show that the kernel of 1 — is 
trivial. Now, it is a standard fact that ker (1 — C u ) = {0} if and only if the associated 
homogeneous RH problem (for say m") has only the trivial solution. But the explicit 
relation between v]/^ spec ^ and m n also establishes a relation between solutions F n and m° 
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Figure 2: Contour for the RH problem for m a . 



of the associated homogeneous RH problems. In view of Proposition 12.31 which states 
that F a = 0, the second step has thus already been accomplished. 

We now establish the above mentioned relation between vj/£ spec ) anc l mQ; derive the 
RH problem satisfied by m a , and finally show that a factorization of v may be chosen so 
that 1 — Cu, is Fredholm with index 0, cf. |25|. 



Let D = {C G C | |C| < 1}. Set 0(C) = |C 3/2 + -< 1/2 and 



m a {() 



*£~W ( c ° -" c c I logC 

1 -i 



for ( 6 % n 



(2.37) 



with {A,-}| =1 being the matrices in Proposition I2.2[ and where n a = 1 if a — 1/2 G No and 

otherwise. By Proposition 12.21 it follows that m a is analytic in D. Let £ = S U <9ID) and 
orient the components of E as in Figure [2J This makes £ a complete contour, meaning that 
C\S can be expressed as the union of two disjoint sets, C\S = f2 + Uf2_, fi + nO_ = 0, such 
that S isjjhe positively oriented boundary of fl + and the negatively oriented boundary of 
Let Ej = %ndD. 

Computations show that m a satisfies the following normalized RH problem. As in 
Proposition 12.21 we use Vj to denote the jump matrix on in the model RH problem. 



Riemann-Hilbert problem for m a 

• m a : C \ X — > C 2x2 is analytic. 

• m Qj+ (C) = m a _(()v(() for ( G E, where 
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A-<7 3 /4J_ ( 1 * 

^ ^ i 1 



C ^C a iogC 



e~ e,T3 t; j e 9,T3j - 







V2 \ - 



A . P 8v3 J_ 
^ 6 V2\-i 



v 0-3/4 



/•- 0-3/4 J_ 



e~ e,J3 A 



-1 



? 1 

m a (C) = I + 0(1/0 asC^oo. 







£0-3/4 
< a logC 



for c e S n D, 

forCe^-nl^e^M}, 
for c e s 3 n D c , 

for(6S 3lJ e{l,3}, 

forCeE i ,je{2,4}. 



The analyticity of m a on S3 H D c follows since #+(C) + $-(C) = for £ < 0. 

It is important to note that v (() — / decays exponentially as ( — > 00 along S. Next 
observe that, at any of the points 0, A, B,C, D of self-intersection of S (see Figure [2]), 
precisely four contours come together. At a fixed point of self-intersection, say P, order 
the contours that meet at P counterclockwise, starting from any contour that is oriented 
outwards from P. Denoting the limiting value of the jump matrices over the jth contour 
at P by v^'(P), we then have the cyclic relation 



„(D(p) ( V (2)(p))-y3)(p) (v^iP))' 1 



I. 



(2.38) 



This is trivial in case P = 0, and follows by direct computation in the other cases. We 
remark that the cyclic relation (12.381) at C is a consequence of the relation (12.121) in the 
case a — 1/2 ^ N , and of (12.131) in the case a - 1/2 e N (see the proof of Proposition 
Q. 

Outside small neighborhoods of the points of self-intersection we choose the trivial 
factorization v + = v, t>_ = / in ( 12.311) . so that oj + = v — I, oj_ = by (12.321) . Using the 
cyclic relations (I2.38p . we are then able to choose a factorization of v in the remaining 
neighborhoods in such a way that uo + is continuous along the boundary of each connected 
component of Q + , and similarly, c<j_ is continuous along the boundary of each connected 
component of 

The exponential decay of v (() — / as ( — > 00 ensures that uo± G £ 2 (S) D L°°(S). From 
this it follows that 1 — is Fredholm in L 2 (S). Indeed, set 



I — vZ , u5+ = v 7_ — I. 
The choice of uj = (£;_,£;+) is motivated by the relations 

UJ_ UJ_ = UJ_ + LU_ , UJ + L0 + = — (i0 + + UJ + ). 

A direct calculation, using C + — C_ = 1 and (I2.40p . shows that 
(1-C U )(1-C S ) = 1 + T, 

where 

Tf = C + ((CL[/(S+ + C_)]) w_) + CL((C + [/(S+ + «_)]) w + ) 



^2.39) 



(2.40) 



f2.41) 



(2.42) 
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for / G £ 2 (£). Standard computations, using continuity of the functions uj + resp. u_ 
along the boundary of each connected component of Q + resp. Q_, show that T is compact 
in L 2 (Y>). Similar computations show that (1 — C%)(1 — CJ) = 1 + 5*, with S compact in 
L 2 (S). So 1 — Czj is a pseudoinverse for 1 — C u , which is therefore Fredholm in L 2 (S). 

It follows from general theory that the index of the operator 1 — C w equals the winding 
number of det v along E, the latter being defined in the natural way. Now, since det v = 1, 
this is trivially zero. This completes the proof of Proposition 12.41 □ 

Remark 2.5 The RH problem for is indeed solvable for all s G C \ P, where T> is 

a discrete set in C (disjoint from R according to Proposition 12.41) . and the solution \|/£f pec ) 
is meromorphic in s with poles in D. To see this, we first observe that the factorization 
(I2.3ip . (12.321) can be done so that u± are both analytic in s. It follows that s i— > 1 — C u is 
an analytic map taking values in the Fredholm operators on L 2 (S). Since we know that 
1 — C u is invertible for s G K, we then get, by a version of the analytic Fredholm theorem 
|44j . that p defined by (I2.34p is meromorphic. Thus m a and hence ty£ pec ' is meromorphic 
in s. 



2.4 Some preliminaries on equilibrium measures 



Before we embark on the steepest descent analysis for the RH problem of Subsection 12. 1[ 
we recall certain properties of equilibrium measures, see [TH HQ]. We use the following 
notation: 

t = ± V t {x) = \v{x). (2.43) 

As explained in the Introduction, we are interested in the case where n/N —>■ 1 as n, iV — > 
oo, which means that we are interested in t close to 1. For every t we consider the energy 
functional /v*(/x) as in (II. 7p . and its minimizer fi t . 

The equilibrium measure du. t = pt dx is characterized by the following Euler-Lagrange 
variational conditions: There is a constant U G M such that 



2 J log | a; — s\p t (s) ds — V t (x) + l t = 0, x G supp p, t , (2.44) 
2 J log | a; - s\p t (s) ds - V t (x) + l t < 0, x G R \ supp p t . (2.45) 

For t — 1, we have that the support of p,y consists of a finite union of disjoint intervals, 
see [13], say 

k 

SUpp/iy = [J[%,6j] 

3=1 

with ai < bi < a 2 < ■ ■ ■ < ak < bk- Due to the assumption that the density pv of p,y is 
regular, we have the following proposition. 

Proposition 2.6 For every t in an interval around 1, we have that the density pt of pt 
is regular, and that supp pt consists of k intervals, say 

k 

supp p t = [J 6j(t)] 

3=1 
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with ai(t) < b\{t) < azif) < • " < dkif) < b^it). In this interval around 1, the functions 
1 1 — > dj(t) and t i — > bj(t) are real analytic with a! At) < and b'At) > 0. 

Proof. See Theorem 1.3 (iii) and Lemma 8.1 of [36J. □ 

For the rest of the proof of Theorem 11.11 we shall assume that supp /iy consists of 
one interval. In the general case (when supp /iy consists of k > 2 intervals) one proceeds 
analogously, but the parametrix away from the end points given in Subsection 12.5.41 must 
then instead be constructed with the help of the ^-function of -B-periods for the two- 
sheeted Riemann surface y 2 = U^ =1 [(z — ctj)(z — bj)] obtained by gluing together two 

copies of the slit plane C \ U^=i [ a j ? bj) m ^ ne standard way (HI [37]. Since the formulas 
will be more complicated in the multi-interval case, but do not contribute to the main 
issue of the present paper, we chose to give the proof in full for the one-interval case only. 

2.5 Steepest descent analysis 
2.5.1 Preliminaries 

We assume from now on that k — 1, so that supp fiy consists of one interval which we 
take as 

supp(/iy) = [a, 0], a < 0. 

Then there is Si > such that fi t is supported on one interval [a t , b t ) for every t G 
(1 — d>i, 1 + d>i), and its density p t is regular. Hence p t is positive on (a t , b t ) and vanishes 
like a square root at the end points, and it takes the form [H] 

Pt( x ) = tt - \l (h - x)(x - a t ) h t (x), for x G [a t ,b t ], (2.46) 

where h t is positive on [at, bt], and analytic in the domain of analyticity of V. In addition, 
h t depends analytically on t G (1 — 5\, 1 + 8\). 

We are going to use the equilibrium measure pt in the first transformation of the RH 
problem. We remark that in [8j [91 [TOj [20] a modified equilibrium measure was used in the 
steepest descent analysis of a RH problem at a critical point. It is likely that we could 
have modified the equilibrium measure in the present situation as well, but the approach 
with the unmodified [i t also works, as we will see, and we chose to use it in this paper. 

In the one-interval case one can show by explicit computation that 

d 4 d b _ 4 

dt 1 t(b t - a t )ht{a t y dt 1 t(b t - a t )h t {b t ) ' 

which indeed shows that -^a t < and ^b t > 0. It follows that b t > for t G (1, 1 + 5i) 
and b t < for t G (1 — 5i, 1). In both cases we have a t < 0. 

We introduce two functions (p t and (p~ t as follows. For z G C \ (— oo, b t ] lying in the 
domain of analyticity of V (which we may restrict to be simply connected, without loss 
of generality), we put 

<Pt(z) = 7i l\{s~bt){s-a t )) l/2 ht{s)ds, (2.48) 
1 ht 



(2.47) 
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and for zG C \ [a t , oo) also in the domain of analyticity of V, 

0t(z) = I !\{s - b t )(s - a t )) 1/2 h t {s) ds. (2.49) 



2 „„ 

It follows from (j23HD that 

¥>t(*) = lV^th t (b t )(z - btf' 2 xt(z\ (2.50) 
where Xt is analytic in a neighborhood of b t and Xt{h) = 1- Taking 



3 



2/3/1 \ 2/3 



/,(*) = ( ^(z)J = l^^j (* - 6*) xT(z), (2-51) 

we see that f t is analytic in a neighborhood of 6 t with ft(b t ) = 0, 
/l \ 2/3 

fl{b t )=l-V=Eiht(b t )J ^0, (2.52) 

and /t(/z) real for real values of z. Hence, in particular, 

/ t (0) > 0, if t < 1, /i(0) = 0, and / t (0) < 0, if t > 1. (2.53) 

Moreover, / t — > /i as i — > 1, uniformly in a neighborhood of 0. We choose a small disc 
£/(°) around and 52 > sufficiently small, so that f t is a conformal map from C/^ onto 
a convex neighborhood of for every t G (1 — 62, 1 + 62). 

Similarly, there exists a disc centered at a < 0, and a £3 > 0, so that 

/3 \ 2/3 

f t (z) = (|ft(*)J (2.54) 

is a conformal map from onto a convex neighborhood of for every t G (1 — 53, 1 + ^3). 

We let Sq = min(5 1 , 5 2 , ^3) and we fix t G (1 — 5q, 1 + 5q). In what follows we also take 
the neighborhoods and U^ a ' as above. 

2.5.2 First transformation Y 1— > T 

We introduce the so-called ^-function: 

g t (z) = J \og(z - s)du. t (s) = J \og(z - s) p t (s)ds, z G C \ (-00, fe t ], (2.55) 

where log denotes the principal branch. Then g t is analytic in C \ (—00, bt]. Define T by 

T(z) = e^ nhai Y(z) e -\ nl ^ e -«J»W^ j zgC\K, (2.56) 

where /t is the constant from (I2.44p -( l2.45p . By a straightforward calculation it then 
follows that T has the following jump matrix vt on R (oriented from left to right): 

( e -n(gt,+ {x)-gt,-(x)) \ x \ 2a e n(5t,+ (a;)+flt,-(a!)-V r «(a;)+Zt)\ 
e n&K,+ (*)-*,-(*)) ) • ( 2 - 5? ) 

Because of the identities, see [HI [TJ] , 

# t)+ (x) + gt-(x) - V t (x) + l t = -2<p t (x), for x > b t , (2.58) 
g t ,+ (x) + g t ,-(x) - V t (x) +l t = -2<p t (x), for x < a t , (2.59) 

we see that the RH problem for T is the following. 
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Figure 3: Opening of a lens around [a t , 0]. 



Riemann-Hilbert problem for T 

• T : C \ R -> C 2x2 is analytic. 

• T + (x) = T-(x) vt(x) for with 



v T (x) = < 



(\ |x| 2a e -2 "-^^^)^ 
1 

e 2n</3 ti+ (x) |a;| 2a 

e 2n<^i,_(a:) j > 

1 



for x < a t , 
for x G (a t ,b t ), 
for x > & t . 



T(z) =I + 0(l/z) as ^ -> oo. 



If n < 0. 1 hen T( z ) = O ( J ) as z -> 0. If a > 0, then T(» = O ( } 1 ) as 2 -> 0. 



2.5.3 Second transformation T 1— > 5 1 

The opening of lenses is based on the following factorization of i>p on (cbti bt) 



v T (x) 



,2nip t ,+ (x) 





\x 



2a 



D 2riLpt,-{x) 



x \-2a e 2nip t ^(x) j 



2a 



— X 



-2a 







1 

I -2a e 2n<pt,+{x) ^ 



Introduce a lens around the segment [a*, 0] as in Figure [3] (recall that a t < 0). In the 
disc £/(°) around we take the lens such that z ^ ( = f t (z) — ft(0), see (I2.5ip . maps the 
parts of the upper and lower lips of the lens that are in into the rays arg£ = 27r/3 
and arg£ = — 27r/3, respectively. Similarly, in the disc LA°) we choose the lens so that 
z \— > C = ft(z), see (I2.54p . maps the parts of the upper and lower lips of the lens that are 
in into the rays arg£ = ir/3, and arg£ = — 7r/3, respectively. The remaining parts 
of the lips of the lens are arbitrary. However, they should be contained in the domain of 
analyticity of V, and we take them so that 

Re (p t (z) < -c < 

for z on the lips of the lens outside and U^ a \ with c > independent of t. 
It is important to note that the lens is around [a t , 0], and not around [a t , b t }. 
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Define S by 





' T{z), 


S(z) = < 


™ (<-. 




, ^ ((-« 










Here the map z i— > (—z) 



1 







1 



for z outside the lens, 

for z in the upper part of the lens, 



for z in the lower part of the lens. 



^2.60) 



(I2.60p and the RH problem for T, we find that S is the unique solution of the following 
RH problem. 

Riemann-Hilbert problem for S 

• S : C \ Ss — > C 2x2 is analytic, where £5 consists of the real line and the upper and 
lower lips of the lens, with orientation as in Figure [3J 

• S+(z) = S-(z) vs(z) for z G Eg, where v$ is given as follows. For t < 1, so that 



&i < 0, we have 



I l z \2a e —2n<f>t{z) 

1 

2a 



\Z 



: 

2a o 



v s (z) = < 







U| 2a e -2n<p t (z) 



I |^|2a g— 2ny>t(z) 

1 

1 

t z \-2a e 2nip t (z) I 



while, for t > 1, so that 6t > 0, we have 

/" 1 |z| 2a e - 2nL Pt(zY 
.0 1 

|^| 2aN 

-|^|- 2q 

3 2rop tl +(z) \z\ 2a 

'l |z| 2Q g- 2?1 Vt(2)\ 

.0 1 J ' 

1 0\ 

/_^\-2q e 2nipt(z) j J 5 



for 2; G (—00, Of), 
for z G (a t ,b t ), 
for z G (6t,0), 
for z G (0, 00), 

for z on both lips of the lens, 



for z G (—00, a t ), 

for 2; G (a t ,0), 

for z G (0,6t), 

for z G (6 t , 00), 

for z on both lips of the lens. 



S'(z) — I + 0(1/ z) as 2 — > 00. 

If a < 0, then S(z) = O ( J jj" ) as z -> 0. If a > 0, then S(z) = O ( } } ) as z -> 



from outside the lens and S'(z) = O 



|z|- 2q 1 
\z\~ 2a 1 



as 2 — > from inside the lens. 
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The next step is to approximate S by a parametrix P, consisting of three parts p(°°\ 
?W, and 

( P(°)(z), for z G EA°) \E 5) 

P( z ) = J pW(z), for 2 G C/^ \Es, (2.61) 

[ PH( 2 ), for zgC \ (£/(«>) U £/(«) U (a t , 0)), 

where £7( a ) and £/(°) are small discs centered at a and 0, respectively, that have been 
introduced before. The parametrices p(°°), p( a ) and P^ are constructed below. 

2.5.4 The parametrix P (oo) 

The parametrix p(°°) is a solution of the following RH problem. 

Riemann-Hilbert problem for p(°°) 

• P(°°) : C \ [at, 0] -> C 2x2 is analytic. 

( ) ( ) ( |x| 2q \ 

x) I , i_2a n ) ^ (fit ? 0) , oriented from left to right. 

• p(°°)(z) = J + 0(l/z) as z -> oo. 

The RH problem for p(°°) can be explicitly solved as in [9]. Take 



-a 



D(z) = z a <p[ — — -) , for z G C\ [a t ,0], (2.62) 



-at 

where (f>(z) = z + (z — 1) 1//2 (z + l) 1 / 2 is the conformal map from C \ [—1,1] onto the 
exterior of the unit circle. Then D + (x) D_(x) = \x\ 2a for x G (a t ,0). It follows that 
P>(oo)~°" 3 P^°°'(z) D(z) aa satisfies the normalized RH problem with jump matrix (_°io) 
on (a t ,0) (oriented from left to right), whose solution is well-known, see e.g. (THE], and 
it leads to 

/ | (J3t(z) + A(z)- 1 ) h{pt(z) - friz)- 1 ) \ 
P (oo) (z) =P>oo) CT3 V fKJ ,lWU ] D(z)~ a:i , (2.63) 

for z G C \ 0], where 

p t (z)=[ — ) , for zgC\ [at,0]. (2.64) 
\z — a t J 

2.5.5 The parametrix P^ 

The parametrix P( a ) is defined in the disc U^ a ' around a, where P^ satisfies the following 
RH problem. 
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Riemann-Hilbert problem for 

. p(a) . jj(a) \ s 5 ^ C 2x2 is analytic. 

• p^\ z ) = P<?\z) v s (z) for z G £/ {a) n Eg. 

• (P^*)) -1 = J + 0(n -1 ), as n -> oo, uniformly for 2 G <9£/( a ) \ Eg. 
We seek p( a ) in the form 

p( >(z) = P»(z) e nMz)a3 (-z)- au \ for z G U {a) \ Eg, 

where (— is defined with a branch cut along [0, 00). Then p(°) satisfies a RH problem 
with constant jumps and can be constructed in terms of the Airy function in a standard 
way; for more details see the presentation in [11] . 

2.5.6 The parametrix P^ 

The parametrix p(°\ defined in the disk around 0, should satisfy the following RH 
problem. 

Riemann-Hilbert problem for p(°) 



• p(0) : [/(o) \ s 5 -»• C 2x2 is continuous and analytic on £/(°) \ Eg. 

• Pj. (z) = P_ (z) vs{z) for z G Eg Pi C/W (with the same orientation as Eg). 

• F(°>(z) (P^ 00 )^))^ 1 = J + O^- 1 / 3 ), asn -f 00, i -> 1 such that n 2 / 3 (t-l) = 0(1), 
uniformly for z G dU^ \ Eg. 

• P(°) has the same behavior near as S has (see the RH problem for S). 

A parametrix p(°) with these properties can be constructed using a solution of the 
model RH problem of Subsection 11.21 The construction is done in three steps. 

Step 1: Transformation to constant jumps. We seek p(°) in the form 



P (o) ^ = p (o) ^ e nM*>s z -«o3 ; for z G C/(°) \ E s , (2.65) 

where as usual z~ a denotes the principal branch. It then follows from the RH problem 
for p(°) that P^ should satisfy the following RH problem. 

Riemann-Hilbert problem for p(°) 



P(°) : C/(°) \ Eg -»• C 2x2 is continuous and analytic on U<® \ E< 
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For z G £<? D o^ -* , we have 



1 ) ' 



pf (z) = P^(z)x{ 



(0), 



1 0\ 

o i\ 
-1 oi ' 



for z G (0,oo) n f/ (0) , 

for z in ?7^°' ) on the upper lip of the lens, 

for z G (-oo,0) nf/W, 



I 

-2am 1 ) > for z in on the lower lip of the lens. 



uniformly for z G dU^ \ £5. 
If a < 0, then 





2 




2 




( 


2; 


a 


2; 


:) 



P (0) (z) =0 
while if a > we have that 
p(°)(z) =0 



as z — > 0, 



P (0) U) =0 






; 


2 






2 






z 


: 






z 



Z| |Z 
z| a 1 2 

— a 

— Q I 



as z — ► from outside the lens, and 
as 2 — ► from inside the lens. 



Note that the jump matrices of do not depend on t. 

The reader may note the similarities between the above RH problem for p(°> and the 
RH problem for ty a from Subsection II. 21 In the next step we show how we can use ty a to 
construct a solution of the RH problem for p(°\ 



Step 2: The construction of P (0) in terms of ty a . Recall that £5 in was taken 
such that z 1— > f t (z) — ft(0) maps £5 PI 

£/(0) 

onto a subset of £, where £ is the contour in 

the RH problem for ty a , see Subsection 11.21 

We choose any solution "$? a of the model RH problem and we define P^ by 

P(°\z) = E(z) ^ a (n%(f t (z) - /*(0));n§/ t (0)) , for z G £5, (2.66) 

where E = E h n is analytic in U^°\ Taking P^ as in (12351) with P(°) as in (I2.66P we find 
that all the conditions of the RH problem for p(°' are satisfied, except for the matching 
condition 



P<°)(z) (P (oo) W) _1 = I + 0(n- 1 / 3 ), 
as n — ► 00, t — ► 1 such that n 2 ^ 3 (t — 1) = 0(1). 



^2.67) 
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Step 3: Matching condition. To be able to satisfy (I2.67P we have to take E in the 
following way 



E(z) = P {oo) (z) z a ° 3 



^ 'n»(/ t (z -/t(0) 



for z G U {0) \ [at,0], 
(2.68) 



where both branches are taken as principal. Clearly then E is analytic in \ [at,0]. 
It turns out that E has analytic continuation to U^°\ This follows by direct calculation, 
but it relies on the fact that we chose [a t , 0] as the jump contour for p(°°). 

With the choice (12.681) for E, we now show that (12.671) is satisfied as well. By (I2.65p . 
§ZM§ , we have for z G dU<® \ E s , 



P®{z) = E{z) ^ a (n 2 Hf t (z) - / t (0));ni/ t (0) ) e n <"°> z 



-OUJ'i 



and we are interested in the behavior as n — > 00, t — > 1 such that n 2 / 3 (i — 1) = 0(1). 
We show first that n 2 / 3 / t (0) remains bounded. 

Lemma 2.7 Suppose n — > 00 7 t -> 1 swc/i tfjcrf n 2 / 3 (t - 1) = 0(1). T/jen n 2 / 3 / t (0) 
remains bounded. More precis ely, ifn 2 /*(t - 1) -> L 6 1, i/ien 



^ 2/3 /*(0) - -c 2 , y L = s, 



where 



C2,v = (ci,y) 



2/3 



d6 t 
dt 



(2.69) 



^2.70) 



t=l 



and c^y zs £/ie constant in ( jl .81) . 
Proof. It follows from f[2T5Tj) . that 



/t(0) 



2/3 



V^fkibt) (-b t )x 2 t /3 (o) 



2/3 



d6 f 
dt 



+ o((t-r 



as i — ► 1. 



By dLHD and (12^61) . we have 



so that fl2T69|) - fl2T70|l indeed follows if n 2 ' z (t — 1) —> L. 



[2.71) 



□ 



If we use the formula (12.471) for the t-derivative of b t at t — 1, then we find from (12.701) 
that 



c 2y = 2(-ay^c-J 3 . 
Now we continue with the proof of (I2.67p . 



^2.72) 
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Lemma 2.8 Suppose that n -> oo ; t -> 1 suc/i n 2/3 (t - 1) = 0(1). T/ien (J2SZD 

Proof. In the proof all O-terms are for n — > oo, t — > 1 such that n 2 ' 3 (t — 1) is bounded. 

By Lemma [2.71 the values n 2 ^ 3 f t (0) remain bounded. Since the asymptotic condition 
(c) in the RH problem for ty a is valid uniformly for s in bounded subsets of R, we find by 

UM, (ESSD, and flXSHD 

P<°>(*) = i?(*) (nHf t (z) - f t (0))y a3/4 ^= (J j) (1 + 0(^/3)) 
x exp (-6(n 2 /\f t (z) - f t (0));n 2 / 3 f t (0))a 3 ) e n ^z~^ 

= pM(z)(I + Oin- 1 ' 3 )) exp (- (d{n 2 ' 3 {f t {z) - / t (0)); n 2/3 /*(0)) - n^) a 3 ) 

(2.73) 

uniformly for z G dU^°\ As before we denote s) = §C 3 ^ 2 + ^C 1 ^ 2 - 

The next step is to evaluate the expression in the exponential factor. We have 

6{n 2 / 3 {f t {z) - f t mn 2/3 W)) - rupt 

= \n [(/,(*) - f t (0)f 2 - (ft(z)f 2 ] + nf t (p)(J t (z) - f t (0))^ 2 . 

We will show that this is 0(n~ l l 3 ) uniformly for z E dU^ . To that end, it is enough to 
show that 

F(t,z) := (f t (z) - f t (0)) 3/2 - (ft(z)) 3/2 + |/*(0) (f t (z) - f t (0)) 1/2 

= 0((t-l) 2 ) as*->l, (2.74) 

uniformly for z e dU^°\ 
By (12.531) . we have 

F(l,z) = 0. (2.75) 
Moreover, 

^Ht, z) = \ (/*(*) - / t (0))i | (f t (z) - / t (0)) - |(/ t (*))* | f t (z) 

+ ~ (j t Mu)) (/*(*)- /*(o))* 

+ ^ / t (o) ~ (/«(*) - / t (o))-s | (/ t (*) - / t ( )). 

Let t — 1 and again use ( 12.531) and (12.511) . Due to cancellations one finds 

^F(M) = 0. (2.76) 

Since, in addition, F(t, z) is analytic in both variables and bounded with respect to z in 
dU^ , it follows from a Taylor expansion that F(t, z) = 0((t - l) 2 ), as claimed in (1277111 . 
Thus 

0(n 2 / 3 (/ t (z) - / t (0)); n 2 / 3 /*(0)) - = 0(n^ 3 ), 
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Figure 4: Contour for the RH problem for R. 

so that (12.73P leads to 

p(°)( z ) = p(p°\ z ) (I + 0(n- 1/3 )) , 

uniformly for z G dU^ . Then (12.671) follows since p(°°Xz) and its inverse are bounded in 
n and t, uniformly for z G dU^°\ □ 

This completes the construction of the parametrix p(°\ 

Remark 2.9 The local parametrix p(°) is constructed with the help of a solution ty a of 
the model RH problem. Since the solution ^> a is not unique (see Proposition 12.11) . the 
local parametrix is not unique. In what follows we can take any p(°) and it will not affect 
the final results (Theorems 11.11 and 11.21) . 

2.5.7 Third transformation S i— > R 

Having p(°°), P<°), and P<®, we take P as in (12.611) . and then we define 

R(z)=S(z)P-\z), for z G C\ («9P (0) U<9P (a) US S ). (2.77) 

Since S 1 and P have the same jump matrices on P (0) n£ s , f/ (ct) nS 5 and (a, 0)\(P (0) UP (a) ), 
we have that R is analytic across these contours. What remains are jumps for R on the 
contour shown in Figure H] with orientation that is also shown in the figure. Then, R 
satisfies the following RH problem. 



Riemann-Hilbert problem for R 

• R : C \ S R -> C 2x2 is analytic. 

• = R-(z) Vr(z) for z G Sfl, where 

r p(°°) (pW)- 1 , ondUW, 
v R =< p(°°) (pW)" 1 , on<9P (a) , (2.78) 
[ w 5 (P^ 00 ))- 1 , on S K \ (9£/(°) U 

• P(z) = 1 + 0(l/z) as ^ -> oo. 

Now let n — > oo, t — * 1 such that n 2//3 (t — 1) = 0(1). It then follows from the 
construction of the parametrices (see in particular the RH problems for p(°) and p( a )) 
that 

I + Oin- 1 / 3 ), ondU®, 
v b = < , (2-79) 
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Furthermore, by regularity of the eigenvalue density, there is a constant c > such that 

Reip t (z) > c> 0, for z G E# PI (0, oo), 

Re (pt{z) > c > 0, for z G £« fl (— oo, a), 

Rep t (z) < -c < 0, for zGE R \ (<9f/ (0) U df/ (a) U R). 

This implies (see the RH problem for S 1 ) that v$ = I + 0(e~ cn ) uniformly on £r\ (dU^ U 
at/( a )), so that by (frZH]) 

ujj = / + 0(e" 2cn ) on S B \ (<9£/ (0) U <9f/ (a) ). (2.80) 

The O-terms in (12.791) and (I2.80p are uniform on the indicated contours. In addition, it 
follows from <|2758|> . fl2"39"|) . (P235D, and tlie growth condition (TO]) on V that for any C > 
there exists r = r(C) > 1 such that <ft(x) > Clogx for x > r, and y5t(x) > Clog |x| for 
x < — r. Combined with (12.801) this implies that 

IK - / llL2( Si;A( o C/ (o) uac7 (a) ) ) = 0(e" 2cn ), as n -> oo. (2.81) 

Thus, by (12?79"|1 - (ICT]1 . as n -> oo and t -»• 1 such that n 2/3 (t - 1) = 0(1), the jump 
matrix for R is close to / in both L 2 and L°° norm on indeed 

IK - I\\L*(x R )nL~(x R ) = 0{n~ 1 ^). (2.82) 

Standard estimates using L 2 -boundedness of the operators C± on L 2 (£#) together with 
the correspondence between RH problems and singular integral equations now imply that 

R( z ) = l + 0(n~ 1/3 ), uniformly for z G C \ E R , (2.83) 

as n — > oo, t — > 1 such that n 2 ^ 3 (t — 1) = 0(1). To get the uniform bound (12.831) up to 
the contour one needs a contour deformation argument. Again, see the presentation in 
jll] for more details. 

This completes the steepest descent analysis of the RH problem for Y. 

2.6 Completion of the proofs of Theorem 11.11 and 11.21 

Having completed the steepest descent analysis we are now ready for the proofs of Theorem 
11.11 and 11.21 We start by rewriting the kernel (12.31) for x,y G fl R according to the 
transformations Y t— > T i— > S h- > R that we did in the steepest descent analysis. To state 
the result it is convenient to introduce B = B n ^ as 

B(z) = R(z)E(z), for z G U {0 \ (2.84) 

where E and R are defined in (12.681) and (I2.77p . We also define for x, s G R the column 
vector 

(M*;s)\ { (J), forx>o, 

^afos) = = < v/ /-,\ (2.85) 

V^W S V # a ,+(x; s)e- Q ™ 3 [ J , for x < 0, 



cf. (11.121) . We then have the following result. 
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Lemma 2.10 Let x,y E U {0) D R. Then, 

K n , N (x,y) = 2m( l_ y) (4(^(/*(2/)-/*(0));nl/ t (0))) T (_° 1 J 

x B-\y)B(x)i/J a (n§(/ t (x) -/ t (0));nf / t (0)) . (2.86) 

Proof. We start from the formula (12. 3p for the eigenvalue correlation kernel. Using 
(I2.56P we obtain, for any x,y G R, 

K nN (x,y) = \x\ a e^^^'+W"^' 1 ^''' M a e i«( 2 9t,+(2/)-Vt(?/)+k) 



x 



2ni(x — 



y^(0 l)T-\y)T + (x)Q. (2.87) 



Using (12.581) and the fact that g t>+ = <? t> _ on (b t , oo), it follows that 2g t — V* + It — —2(pt 
on (fe t , oo). Then, by analytic continuation, 2g tj+ — Vt + k = —2(p t ,+ on all of R. Therefore 
we can rewrite (I2.87P as 

K n , N (x,y) = \x\ a e- n ^ | y |« c ^ + W_L_- (0 l) T+ l (y) T + (x) Q . 

(2.88) 

Now we analyze the effect of the transformations T t— > S* i— > R on the expression 
| x |a e -^t,+(s)T + (x) (J) in case x G ?7 (0) H R. The result is that for x G U {0) n R, 



m e 



*' +W T + W (J) = B(aO*«,+ (ni(/ t (ar) - / t (0)); n§ / t (0)) Q (2.89) 



in case a; > 0, and 



\x\ a e- n ^+^T + (x) Q = B(z)* a , + (n§(/ t (a;) - / t (0));n*/ t (0)) e^™ 3 Q 



(2.90) 



in case x < 0. Since the calculations for (12.891) are easier, we will only show how to obtain 
). If x G £/ (0) n R and x < 0, then it follows from fl230"|) that 



V^ (l) T + (,) = \x\ a e- n ^S + (x) ( lx , 2a ^ +(xi 



= S + (x) (|x| Q e- n ^- +(:c) ) <T3 (J\ . (2.91) 

From (EZZD, (123TD . (12351) . (12381) . and (l2~84l) . we find that 

S+(x) = B(x)* a>+ (ni(f t (x) - / t (0));ni/ t (0)) (e^x^ . 



Inserting this into (I2.9ip and noting that x + a |x| a = e am we indeed obtain (12.901) . 
In a similar way, we find for y G n R, 

lyl a e -nv tM v) (o i) T -\y) = (0 l) ^ ( n i(/ t (y) - /*(0));n§/ t (0)) iT^y), 

(2.92) 
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in case y > 0, and 

\ yre -n Vt Mv) (o 1) T -\y) = (-1 1) ( n l(/ t (y) - / t (0));n§/ t (0)) fl-fy), 

(2.93) 

in case y < 0. To rewrite (I2.92p and (I2.93P we use the following fact, which is easy to 
check. If A is an invertible 2x2 matrix having determinant 1, then 

If we apply $ZM) to in fl2T92|) and fl2T93|) . then we get 

lyl a e -n VtM v) (o i) T -\y) = (1 0) (n§(.f t (y) - / t (0));n§/ t (0)) 

LT^y), (2.95) 



1 
-1 



in case y > 0, and 
1^ 



e -n^ + (y) (o i) r- 1 ^) = (1 1) e- Q ™ 3 < + ( n »(/ t (y) - /*(())); n*/^ 

B~\y), (2.96) 



1 
-1 

in case y < 0. Then fl2~86l) follows if we insert fl2T89|) . fl2T90|) . fl2T95|) . and fl2T96|) into (12T88D 
and use the definition (I2.85p . □ 

As in Theorem 11.11 we now fix x, y G R. We define 

Xn= (^)^' yn= (^^ (2 ' 97) 

where Ci is the constant from (11.81) . 

In order to take the limit of (cin)~ 2 / 3 K nt N(x n ,y n ) we need one more lemma. Recall 
that B = RE is defined in (l2~Mj) . 

Lemma 2.11 Lei n — > oo, t — > 1 st/c/i i/iat n 2 / 3 (t — 1) — > L. Let x, y G M an<i Ze£ x n and 
y n defined as in ( I2.97p . TYien £/ie following hold. 

(a) n 2 /3 /t(0 ) - S; 

(b) n 2 / 3 (/ t (x„) - / t (0)) - x and n 2 l\f t {y n ) - f t (0)) - y, 

(c) £? _1 (y n ) B(x n ) =1 + where the implied constant in the O-term is uniform 
with respect to x and y. 



Proof, (a) This follows from Lemma [2.71 

(b) By dUBD and (12^461) we have c x = \^f=ahi(ti), so that f[(0) = c/ 3 by fl232|) . 
Taking note of the definitions (I2.97p . we then obtain part (b), since ft — > /i uniformly in 
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(c) We have 
R-\y n ) R(x n ) = 1 + R-\y n ){R{x n ) - R{y n )) 



I + R-\y n )(x n -y n ) [ R'(tx n + (l-t)y n )dt. (2.98) 

Jo 



Recall that R is analytic in f/ (0) , and that R(z) = I + 0(n~ 1/3 ) by ( |2~83j) . uniformly in 
Since det R = 1, we find that R 1 (y n ) remains bounded as n -> oo. It also follows 
from (I2.83P and Cauchy's theorem, that R'(z) = 0{n~ 1 ^) for z in a neighborhood of the 
origin. By (12.981) we then obtain 



R-\y n ) R(x n ) = I + Ol X —^- ) . (2.99) 



n 



Using analyticity of E in a neighborhood of the origin with E(z) = Oin 1 ^), see fl2T68|) 
and the fact that det E = 1, we obtain in the same way 



E~ l (y n ) E{x n ) = 1 + [^rjrj • (2.100) 

The implied constants in f!2.99|) and (I2.100p are independent of x and y. 

Using fl2T99|) . fl2TT00|) . and the fact that £(ac n ) = O^ 1 / 6 ) and E- l (y n ) = 0{n l / & ), we 
obtain from ( 12784]) 

5- 1 (y n ) J B(x n )= J E- 1 (y n ) fl + o(^y\ E(x n ) 

= E-\y n ) E(x n ) + 0(n 1/6 ) O (j^j <V /6 ) 



I + O 



x-y 



n l/3 

This completes the proof of part (c). □ 

Proof of Theorems 11.11 and II. 2L We let n, N — * oo, t = n/N — > 1, in such a way that 
n 2 / 3 (t — l) — > L. Then by parts (a) and (b) of Lemma [2.111 we have 

j a (n 2/3 {ftM - / t (0)); n 2 / 3 f t (0)) -> ^(x; s) 

and similarly if we replace x n by ?/ n . The existence of the limit (11.91) then follows easily 
from Lemma 12.101 and part (c) of Lemma 12.111 which proves Theorem 11.11 
We also find that the limiting kernel K^f 9e (x,y; s) is given by 

K °"^ « s) = mhv^ a)T (-°i J) te s) 



and so f 1 1 . 1 3 j) follows because of (I2.85p . The model RH problem is solvable for every s £ K 
by Proposition 12.41 and so we have also proved Theorem 11.21 □ 



CRITICAL EDGE BEHAVIOR IN RANDOM MATRIX ENSEMBLES 



33 




FN 



FN 



Figure 5: Contour for the RH problem for \l/ 



FN 
v 



3 Proof of Theorems 11.41 and 11.5 



We prove Theorem 11.41 and Theorem 11.51 by first establishing, with the help of [2], a 
connection between the model RH problem and the RH problem for Painleve II in the 
form due to Flaschka and Newell [22J. We can then use known properties of the RH 
problem for Painleve II to prove the theorems. 



3.1 The Painleve II RH problem 

We review the RH problem for the Painleve II equation q"(s) = sq + 2q 3 — v, as first given 
by Flaschka and Newell [22], see also [21] and [25]. We will assume that 

v > -1/2. 

The RH problem involves three complex constants ai, a^-, a$ satisfying 

cli + a2 + a3 + Gliomas = — 2i sin utt, (3.1) 

and certain connection matrices Ej. 

Let Sj = {w G C I ^vr < argu; < ^±tt} for j = 1, . . . , 6, and let £™ = C \ [J ■ Sj. 
Then T, FN consists of six rays Tif N for j = 1, . . . , 6, all chosen oriented towards infinity 
as in Figure El The RH problem is the following. 



Riemann-Hilbert problem for 

• V™ : C \ Z FN -> C 2x2 is analytic, 

*S? = (J on , 



34 



A.R. ITS, A.B.J. KUIJLAARS, and J. OSTENSSON 



= <- ( J T ) on 



• *f ^(u;) = (I + 0{l/w))e- i ^ w!i+sw ^ as w -> oo. 

• If 1/ - | ^ N , then 

*fW = %)(^ jl)^', fortoefi,-, (3.2) 
where I? is analytic. If z/ G | + No, then there exists a constant k such that 

(w) = #0) o ) % for w G 5 is (3.3) 

where B is analytic. 

The connection matrix E\ is given explicitly in [2H Chapter 5]. It is determined (up 
to inessential left diagonal or upper triangular factors) by v and the Stokes multipliers 
ai, a 2 , and a 3 , except in the special case 

u = ^ + n, a t = a 2 = a 3 = i(-l) n+1 , n G Z, (3.4) 

where an additional parameter c G CU {oo} is needed. For example, for v G" \ + No and 
1 + a±a 2 7^ 0, we have 



'd \ / 1 -r— 

1 + aia 2 

1 + a\a 2 e vm + a 2 

V 2 COS Z/7T 2 COS Z/7T / 



v rf 



-i 



(3.5) 



where d ^ is arbitrary. In the special case 03.4)) . when £7 X depends on the additional 
parameter c G C U {oo}, by [2U Chapter 5, (5.0.21)] we may take E\ as 

E x = r f\, if c G C, while Ex = K if c = oo. (3.6) 

Assuming that the branch cuts for the functions in (13.21) and (13. 3 j) are chosen along 
argw = — 7r/6, we obtain the other connection matrices from Ex through the formula 

E j+1 = Ejvf N , j = l,...,5, (3.7) 

where vJ N is the jump matrix on ■ We shall refer to the Stokes multipliers a 1; a 2 , 
and 03, and in the special case 03.4)) also to the additional parameter c, as the monodromy 
data for Painleve II. We note that in the special case (13.41) we have k = in 03.31) . 
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The special case ( 13.41) has geometric interpretation. Indeed, (13.41) describes the singular 
point of the algebraic variety (13.11) . that is, the point at which the (complex) gradient of 
the left-hand side of ( 13.11) vanishes. The singularity may be removed by attaching a copy 
of the Riemann sphere (see also [30]). 

The monodromy data does not depend on s. The RH problem is uniquely solvable, 
except for a discrete set of s- values, and its solution d e p enc i s on s through the asymp- 
totic condition at infinity. We write ty^ N (w; s) if we want to emphasize its dependence 
on s. If we take 

q(s) = 2% lim w (V™(w; s)) 10 e i ^ wi+sw)a \ (3.8) 

then q satisfies the Painleve II equation q" = sq + 2q 3 — v. In addition ^f^ N satisfies the 
Lax pair for Painleve II 

a f-Aiw 2 - i(s + 2q 2 ) 4wq + 2ir + \ 

JT* = LV, L=[ ), (3.9) 

dw V Awq-2ir + ^ Aiw 2 + i(s + 2q 2 )J 

l*=p% p=(~: w i), (3.io) 



ds \ q iw / 

where q = q(s) and r = r(s) = q'{s). In this way there is a one-to-one correspondence 
between monodromy data and solutions of Painleve II. 
We also need the more precise asymptotic behavior 

S » - {' + h Cw -Hi)) + °< 1 /- 2 )) e" 1 ^' («D 
as w — > oo, where 

H(s) = (q'(s)) 2 - sq\s) - q\s) + 2vq{s) (3.12) 

is the Hamiltonian for Painleve II. Note that H' = —q 2 . 

We finally note that ty„ satisfies the symmetry property 



V 

^ N (w; s) = ^ N (-w; s)^, (3.13) 

where a\ = (?q). Indeed, by a straightforward calculation (see also [2H Chapter 5]) 
we check that the function ai^^ N (— w; s)ai solves exactly the same RH problem as the 
function ty„ N (w] s). Unique solvability of the RH problem yields equation (13.131) . 



3.2 Connection with ^ 



a 



The Hastings-McLeod solution of Painleve II corresponds to the Stokes multipliers at = 
—e um , ci2 = 0, and 03 = e~ un \ This is not the solution that interests us here. We use 
instead the solution corresponding to 

a 1 = e- u ™, a 2 = -i, a 3 = -e uni . (3.14) 

For these Stokes multiplies (I3.14p we obtain from (13.51) the following connection matrix 
Ei in case v | + N (where we take d = (e um — i)/{2 cos uix)) 

e v ™-i ie uwi + V 

Ei = J 2COSZ/7T 2COSZ/7T I . (3.15) 

-e~ v ™ 1 
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For v G | + 2N , it follows from 03.141) and the formulas in [Ml Chapter 5, (5.0.18)] that 
we can take 

* = ?) • (3-16) 

If i/ G | + 2N , then we are in the special case (13.41) . We then choose c = i, so that for 
v E \ + 2No we have monodromy data 

fll = e"^ = -i, a 2 = -i, a 3 = -e vm = -i, c = i. (3.17) 



Lemma 3.1 For any v > —1/2, we have that 

(E 2 ) 21 = (E 3 ) 21 = 0. (3.18) 



Proof. In all cases we may check from (13.141) . (I3.15p . (I3.16p . (13.171) . and (13. 6p that the 
second row of E x is given by (— a\ l). So by (13. 7p we have that E 2 — Ei ( \ ) is upper 
triangular. Then also E 3 = E 2 (J a £ ) is upper triangular and therefore (I3.18P holds. □ 



The following proposition holds for more general monodromy data, and it was estab- 
lished in [2], see also [35]. For the reader's convenience we present a detailed proof for 
our particular case. 



Proposition 3.2 ([2]) For a > —1/2, let ^ 2 a+i/2 ^ e un ^Q ue solution of the RH prob- 
lem for Painleve II with parameter v = 2a + 1/2 and monodromy data (13.141) in case 
a ^ No (so that v ^ | + 2No), and monodromy data (13.171) in the special case a £ No. 
Then, for any r\ = r](s), we have that 

* Q (C; s) = (J 8) J) C- CT3/4 -^ (J {) e™* ^ a N +1/2 (w; -2 1 /3 s)e — 3/4 (3 19) 

where w = e 7r */ 2 2 -1 / 3 ^ 1 / 2 with Imw > 0, is a solution of the model RH problem for 
given in Subsection \l. 6 A 



Proof. Because of Proposition (12. ip we may take n(s) = without loss of generality. 

Clearly is analytic on C\E. The correct asymptotics as £ —>■ oo follows immediately, 
as well as the correct jumps across Si, £2, and £4. A little bit more work is needed to 
check the jump across S3 = (—00, 0) and the behavior at z = 0. 

In order to analyze the jump across S 3 , we suppose that ( G S 3 . Then we have that 



= e ^ 2 2^ 3 Cl /2 = -w- = -e" /2 2- 1/3 ( V2 (< 0), 
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and hence by (I3.19P and the symmetry property (13.131) . 

WC; s) = CT 74 ^ (J j) e-/X: +1/2 (%; -2 1 /3 s)e — 3/4 

=c ct3/4 (-; °)^() ;)e^3/4^ r Ar i/2( _ w _._ 2 i/3 S)e -^3/4 

_ ^-^3/4 J_ / 1 \ 7ri CT3 /4^FAr / . _ ? l/3 \ ( M 

-U ^2^j W 2a+1/2 ^_, ; sje ^ 

_ ^-a 3 /4 J_ / 1 A iCT3/4 FA r , _ _ i /3 , ^3/4 (0 A 

-U ^2\i \) e v 2a+iM«>-> l s)e Q j 



*a,-(C;s) 



1 
-1 



B(w) 



This shows that has the correct jump across S3, and it follows that ty a satisfies the 
parts (a), (b), and (c) of the model RH problem. 

Consider now a neighborhood of the point ( = 0. We recall that (13. 2p or (13.31) holds 
with B(w) analytic at 0. A corollary of the symmetry property (I3.13P is the equation 

if v £ \ + No, 
as w — > 0, if v e \ + N , 

which yields the formula (cf. [24], Chapter 5]) 
B(0)a 3 = a x B(0). 

The last relation, together with det B(0) = 1, in turn implies that -B(O) can be represented 
in the form 




If C e flj then w E S n ^), j = 1, 2, 3, 4, where n denotes the permutation 



7T 



12 3 4 
3 4 12 



Therefore, for the function ty a ((] s) defined by equation (13.191) with rj(s) = 0, we find that 
(assuming that a ^ ^Nq) 



tf a (£; S ) = C~ CT3/4 ^| f • J J e ™3/4 S ^) e -™ 3 /4 (j + 0(C 1/2 )) C CT3/4 C a<T3 ^7r 

= rCT3/4 K! T) (0 6 °x)( J +°(c 1/2 ))r 3/4 r-K 



C -s/4 (J + 0(C l/ 2) ) r 3/4 C a. 3 A M ^ 

0(l)C a:i (l asC-Oin^, (3.20) 
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where we have introduced the notation 

, (2a+l/2)<7 3 



Ej = e m ^ (e"^" 1 / 3 ) E je -^\ (3.21) 

From (13.201) it immediately follows that ^(C; s) = 0(C~' a ') as ( — > 0, which is the 
required behavior in the model RH problem if a < 0, or if a > and j G {2, 3}. If a > 
and j G {1,4}, then n(j) G {2,3}, and it follows from Lemma I3TT1 and (I3.2ip that 



Then (13.201) also yields the required behavior of ^ a ((; s) as ( — * in Qi U fi 4 . 

The calculation leading to (I3.20p is valid for v ^ | + No, or a ^ |No- In fact it is also 
valid if a G No, since then we are in the special case (13. 4p where k = in (13.31) and so 
no logarithmic terms appear. Logarithmic terms only appear if a G ~ + No, and then a 
similar calculation leads to 



with £y again given by (13.211) . Since a > 0, the required behavior as £ — > then follows 
in a similar way. 

This completes the proof of the proposition. □ 
3.3 Differential equation 

Recall that *ff„ has the Lax pair (13.91) and (I3.10p . Then defined by (I3.19P also satisfies 
a system of differential equations. It will involve the solution q of the Painleve II equation 
with parameter v = 2a + 1/2 and monodromy data (13. 14j) or (13.171) . We put r = q' and 

U(s)=q 2 (s)+r(s) + ^ (3.22) 
V(s) = q 2 (s)-r(s) + S -. (3.23) 

The functions U and V both satisfy the Painleve XXXIV equation in a form similar to 
dTTSD , namely (cf. [22 Chapter 5]): 

^M = W + M ™-"M-S- (3 ' 24) 
v.(.,-^ + ^.,-.K(.,-^£. 

Then we obtain the following differential equations for 
Lemma 3.3 Le£ ty a be given by (13.191) . 

(a) lfr] = 0, then satisfies 

iU a (C;s)=A* a (C;s), (3.26) 
JU a (C;s)=fl*«(C;s), (3.27) 
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where 



^=1 .... Zu*J n .L m,. J,, J'l. ( 3 - 28 ) 



-< + «2" 1 /3y(_2 1 /3 s ) 2 1 /3 g (_ 2 V3 



c 




. 2 V3 g (_2 1 /3 s ) 

5=1 , ; I . (3.29) 

(b) For general rj we have that ty a satisfies 

=(*.) ?W-,M i)*° (C;s) ' (330) 

tuzi/i A given fry (13. 28[) . 

Proof. This follows by straightforward calculations from (I3.9p . (13.101) . and (13.191) . □ 

The Lax pair (I3.26p - (I3.27I) . after the replacement £ i— > £ — s, becomes the Lax pair 
from [2J [35] . Equations (I3.24I) - (I3.25I) can also be derived directly from the compatibility 
conditions of the Lax pair f!3.26j) — (I3.27j) in a usual way. 

It is a fact |33J, that the solution q of the Painleve II equation (with parameter 
v = 2a + 1/2 and monodromy data (13.141) or (I3.17P ) has an infinite number of poles on the 
positive real line, see also (I4.29P below. If — 2 1 / 3 s is such a pole then ^20+1/2 ('> ~~ 2 1 / 3 s) 
does not exist. So to be precise, if we assume that rj is analytic on R, then (13.191) does not 
define ty a for values of s G M. which belong to the discrete set of values s where q(— 2 1 / 3 s) 
has poles. 

The relation (13.191) defines ty a for all s e M only if we are able to choose rj so that all 
the poles on the real line of the right-hand side of (I3.19P cancel out. Such a choice of 77 
would require r] itself to have poles at the poles of q(— 2 1 / 3 s). 

We will describe two special choices for rj. The first choice is such that (I3.19P is equal 
to the special solution fya Pec \ which is characterized by the asymptotic condition (I2.14p . 
From Proposition 12.41 we know that \^^ pec ^ exists for all s G R, so that we can already 
conclude that the special choice r] = r/( s P ec ) will have poles at the poles of q(— 2 1 / 3 s), and 
that the real poles of the right-hand side of (13.191) will indeed cancel out. 

The second choice of rj is made so that the differential equation (13.301) takes a nice 
form. It will lead to the differential equation (ll.2ip for ipi and ip 2 - This rj is denoted 770, 
and it is defined by the simple formula 

r !0 ( S ) = t2 1 / 3 q(-2 1 / 3 s), (3.31) 

from which it is already clear that it has poles at the poles of q(— 2 1 / 3 s). For the choice 
(I3.3ip we can already check that the differential equation (13.301) leads to 

^* a (C;s)=4o*a(C;s) (3-32) 

where 

A„=( l "\A fl ° 



V \-Vo 1 



(a + iurj )/C i-iu/C \ 

-< + i(v + ri) +ri (2a + iu7] )/C -(a + iurj )/CJ ' 
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and 

u(s) = 2- 1/3 U(-2 1/3 s), (3.34) 
v(s) = 2- 1/3 V(-2 1/3 s). (3.35) 

3.4 Special choice r]( spec ^ 

Lemma 3.4 Let H be the Hamiltonian for Painleve II as in (13. 12[) . with parameter v = 
2a + 1/2, and let 

v ^ ec \s) = i2~ 2 ' 3 (g(-2 1 / 3 s) + H(-2 l / 3 s)) . (3.36) 

Then the choice rj = r]( spec ^ in (13.191) leads to the special solution \fr^ pec ^ of the model RH 
problem characterized by (12.141) . 

Proof. It follows from ( 13. Ill) and ( 13. 19p by straightforward computation, that 
*«(C; s) = (J fj C-- /4 ^ (J j) e-^ a N +1/2 (w; -2 1 /3 s ) e —/4 



1 
r/ 1 



1 0\ 1/0 i2~ 2 l z {H-q){-2 1 ^s] 



+ C^ 4 ^ Q j) e^/ 4 0(l/C)e- 7riCT3/4 ^= C 3/4 



x C^ 4 ^ (] JJ e -(IC 3/2 +< 1/2 )- (3.37) 

as C — > oo. From (I3.37j) it is clear that we need to take rj = r]^ spec?l in order to be able to 
obtain (12. 14[) . Thus the lemma follows. □ 

From the calculation (13.371) we also note that for any solution ty a of the model RH 
problem we have 

,§(c^ +s cv >3 J_ ( 1 -A r»A __ ^ 2/3 (H-g)(-2^ s ) , 



;jc 3/4 J = - — v ^ ' +o(r i/2 ) 

(3.38) 

as ( — > oo. This property will be used later in the proof of Theorem 11.41 

Since the left-hand side of (13.381) is analytic in s for s 6 R, it also follows from (13.381) 
that H — q does not have poles on the real line. This and similar properties are collected 
in the following lemma. Recall that U is given by (13.221) . 

Lemma 3.5 The following hold. 

(a) H — q has no poles on the real line. 

(b) U has no poles on the real line. 

(c) U has a zero at each of the real poles of q and Uq has no poles on the real line. 

(d) Uq takes the value v — 1/2 at each of the real poles of q. 
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Proof, (a) We noted already that part (a) follows from (13.381) . 

(b) Since H' = — q 2 , we have that 

U(s) = q 2 (s) + q'(s) + s/2 = -(H - q)'(s) + s/2, (3.39) 

and so it follows from part (a) that U has no poles on the real line either. 

(c) Differentiating (13.221) . we obtain 

U' = 2qq' + q" + \ = 2qq' + sq + 2q 3 - v + \ = 2Uq - v + \. (3.40) 

Thus also Uq has no poles on the real line, which means that U has a zero at each of the 
real poles of q. 

(d) Using (13^01) . we get 

(Uq-u+ ~)q = (U' - Uq)q = (Uq)' - U(q 2 + q') = (Uq)' - U(U - s/2). (3.41) 

Since the right-hand side of (I3.4ip is analytic on the real line by parts (b) and (c), we 
conclude that Uq — v + \ has a zero at each of the real poles of q. This proves part (d). □ 

It is well-known and easy to check that each pole of q is simple and has residue +1 or 
— 1. Indeed, the Laurent series for q at a pole s has the form 

q(s) = q ^ + q x (s -s ) H , 

s- s 

where g_i G { — 1, 1}. Using this, one easily verifies that either q 2 + q' or q 2 — q is analytic 
at s (depending on the sign of the residue q~i). Our result that U = q 2 + q' + s/2 is 
analytic on M can then also be stated as follows. 

Corollary 3.6 The solution q of the Painleve II equation with parameter v = 2a + 1/2 
and monodromy data (13.141) or (13.17!) has only simple poles on the real line, with residue 
+ 1. 



3.5 Special choice t]q 

As already announced we will also use the special choice rj = rj Q given by (13.311) . 
By fl3~3TD and <^M> we have that 

Vo(s) - V (spec) (s) = i2- 2 ^ {q(-2^s) - H^'h)) , 

and so it follows from part (a) of Lemma [3.51 that 770 — r/ spec ) is analytic on the real line. 
Since \l/o Pec ^ exists for all s 6 R, it follows that the solution of the model RH problem 
associated with i]q exists for all s G M as well, and it is analytic in s. 

The differential equation for fy a with 77 = 770 is given by (13.321) with Aq as in (I3.33[) . 
It then follows that Aq is analytic on the real line, and we will explicitly verify this by 
rewriting its entries in terms of the function u from (13.341) 



u(s) = 2- 1 / 3 f/(-2 1 / 3 s). 
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The analyticity of u is immediate from (I3.34p and part (b) of Lemma [3751 The analyticity 
of ur] follows from (I3.34p . (13.311) and part (c) of Lemma [3751 Using also (13.401) we get 

u' = 2iur]o + u -l/2 = 2iur]o + 2a. (3.42) 

Next, it follows from (j3722]) . (|3723j) . 3333), (137351) . and (137311 that 

u(s) +r/ (s) 2 = -u(s) -s. (3.43) 

We can use (13.421) and (13.431) to eliminate rjQ and v from the entries in Aq, and we get 
from (137331) that 

\-iC ~ <(« + a) - *(K) 2 - (2«) 2 )/(4<) -u'/(2C) J ' 1 J 

3.6 Proof of Theorem 174] and [T751 

After these preparations the proofs of Theorems 11.41 and 11.51 are short. 

Proof of Theorem 11.41 From (I3.24p and (13.341) it follows that u satisfies the Painleve 
XXXIV equation in the form (1X71^ . 
From (137381) it follows that 



lim 



12J 



which in view of ( I3.39P and (I3.34p leads to (11.201) . This proves Theorem 11.41 □ 

Proof of Theorem 11.51 Let ^> a be the solution of the model RH problem given by 
(I3.19P with rj = r] as in (13.311) . Then 

^* a (C;s)=A)*a(C;s), (3.45) 

with Aq given by ( 13.441) . The differential equation ( I3.45P is valid for ( 6 C \ S. We can 
take the limit ( — > x with x e M\{0} to obtain a differential equation for ^/ a + (x; s), with 
the same matrix A (but with ( replaced by x). Using (I1.12p . we obtain the differential 
equation (11.211) for ipi and ip2- This completes the proof of Theorem 11.51 □ 



4 Concluding remarks 

4.1 The case a = 

The case a = is classical and well understood. We know that KQ d9e (x,y; s) is the 
(shifted) Airy kernel, see (11.111) . We will show here how this follows from the calculations 
from the previous section. 

In the special case a = 0, we have v = 1/2, and then the Painleve II equation has 
special solutions built out of Airy functions. To be precise if Ai and Bi are the standard 
Airy functions, then for any C\ and C2, not both zero, we have that 

q( s ) = £ log (C 1 Ai(-2- 1 / 3 S ) + C 2 Bi(-2- 1 / 3 S )) (4.1) 
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is a solution of q" = sq + 2q 3 — |. These are exactly the solutions that correspond to the 
special Stokes multipliers cti = a 2 = a 3 = —i. The corresponding solutions to the RH 
problem were given by Flaschka and Newell [22j Section 3F(iv)]. For example, for w in 
sector Si we have (see also [24j Chapter 11]) 

^n (w . s) _^o w i/2 (l-iq(s)/w -2-VH/w\(Ai(z) Bi(z)\(-i l\ 

with z = —2 2 ^w 2 — 2 -1 / 3 s and a Q = 2 1 ^^e ln ^ 4 . The expressions for \Pf N (w; s) in the 
other sectors follow by multiplying (14. 2 p by the appropriate jump matrices. 

It follows from (14. ip and (14.21) that the extra parameter c in the monodromy data for 
flUD is 

^ - C 2 

So if we take c = i as in (13.171) then C 2 = and the corresponding solution (14.11) is 

?(s) ^ log Ai(-2-V3 sH _ 2 -V3g^. (4 . 4) 

Note that the solution (14.41) is special among all solutions (14.11) in its behavior for s — > — oo. 
Indeed, from the asymptotic behavior for the Airy functions it follows that for (14.41) we 
have 

q(s) ~ -v / 2(— s) 1//2 as s ^ — oo, 

while for the other solutions (14. ip we have 

q(s) ~ — -v / 2(— s) 1//2 as s — >• — oo. 

So according to Proposition 13.21 we should be using q given by (14.41) and then define 
\l/o as in (I3.19p . If ( is in sector il^, then w = e m ' 2 2~ l ' z C > 1 ' 2 is in sector Si, so that by 

^^ 2[W ' 1 s) ~ 2 c vi-^^r 1 / 2 c" 1/2 ) V Ai '(c + s ) Bi '(c + s)j v 1 -v 

where r? (s) = z2 1 / 3 g(-2 1 / 3 ) as in (T373T]) . Then fl3TT9|) with 77 = r? yields for C G f2 3 , 



e 7TiCT 3 /4 



i+^ ( S )r 1/2 -r 1/2 V Ai ^ + s ) Bi(c + S )\A z 
1 - i V o(s)c- 1/2 r 1/2 A Ai '(c + ^) Bi'(c + S )J^ 1 

W4 /Ai(C + s) Bi(C+s)\ A i 
^2 e VAi'(C + s) Bi'(C+ S )j ^ 1 
A ( Ai(C + s) + zBi(C + s) -(Ai(C + s)-iBi(C + s)) 
72 V-*(Ai'(C + s ) + iBi '(C + *)) i(Ai'(C + s) - zBi'(C + s)) 



— nia 3 /4 
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Since (see e.g. formula (10.4.9) in [lj) 

Ai(z) ± iBi(z) = 2e ± ™ /3 Ai(e T2 ™ /3 z) 

we can write \l/o m the more familiar form 

_ /- / e -/3 Ai (e-^/3 (c + s)) - e -W3Ai( e «/3 (( + s)) \ 
w ^, sj - V2vr ^_ ie -W3Ai / ( e -^/3(C + s)) 2e™/ 3 Ai'(e 2 ™/ 3 (C + s)) y ' ; € 3 ' 

(4.5) 

For ( G Hi we find in a similar way (or by multiplying (14.51) on the right by ( J)), that 

* (<r- s i - T Ai(c + s) ^/3 Ai ( e -^/3 (c + s)) \ 

*o(C,^-V27r^_. A . /(c + s) _ ie -W3Ai'( e -2-/3( C + s ))J> forCe^. (4.6) 

Then it follows from (11.121) and (14. 6 p that for x > 0, 

V>i(x; s) = \/27rAi(a; + s), ^(z; s) = -V^7rzAi'(a; + s), (4.7) 
and a similar calculation shows that (14. 7p also holds for x < 0. Therefore, by (11.131) 

Ai(x + s)Ai'(y + s) - Ai'(x + s)Ai(y + s) 
x - y 

which is indeed the (shifted) Airy kernel. 
4.2 The case a = 1 

The case a = 1 can be solved explicitly in terms of Airy functions as well. Let \l/o be a 
solution of the model RH problem with parameter a = 0. Then for any matrix X = X(s), 
it is easy to check that 

* x (C;s) = (/-iA;( a ))*o(C;a) (4-9) 

satisfies the conditions (a), (b), and (c) of the model RH problem for a = 1. For a special 
choice of X we will have that the condition (d) is also satisfied. 

Let's take \l/o given by (14.61) for ( G fij. Then the condition (d) of the model RH 
problem yields the following condition on X 

(I - \x(s)) ( fl^fl,) = 0(C), as C - 0. (4.10) 



(4- 



( w; \-iAi'(( + s 
The condition (14.101) is satisfied if and only if we take 

Note that the denominator in (14.111) cannot be zero for s El. Indeed, its derivative is 
— Ai(s) 2 , so that it is decreasing for sGl, and since the limit for s — > +oo is equal to 
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0, it follows that Ai'(s) — sAi(s) 2 > for all sGR. Note also that if we take the limit 
x, y — ► in (14. 8p . then 

A^ e (0, 0; s) = Ai'(s) 2 - sAi(s) 2 . (4.12) 

Using (Q2]| . (jH) , (Oil . pi| . and (Q2]l . we obtain that 

, / s / — . / >. / — Ai(x + s)Ai'(s) - Ai(s)Ai'(x + s) , , . 

x; s) = v^Ai (x + s - V2^— /. ' t^/t^ -Ai(s) 

v ' y ' x(Ai'(s) 2 - sAi(s) 2 ) v ' 



K d9e (x,0;s) 
< 9e (0,0; S ) 



^ I Ai(x + S ) - ,i ge ::\ ; Ai(3) 



, . , , — . .. . , — Ai(x + s)Ai'(s) — Ai(s)Ai'(x + s) . ., . 
M x; s) = -^{x + s) + ^ x [ M ^_ a ^) Al (g) 

- y ™, I Ai'(x + s)- K \ e ( X, °' ,s) Ai'(s) 
V ' K e d9e (0, 0; s) 



Thus 

Kr ae (x,y; s 



edge, , ^ifo S^lG/J s) - 1pl(x; 8)lfa(.V, *) 



2ni{x — y) 



-edge, ,,. ,, Kf ge (x, 0; s)Kl dge {y, 0; s) 



To compute the relevant solution w of the Painleve XXXIV equation for a = 1, we 
may assume that we have taken ^/^ pec in (14.91) . and then use (ll.20p . (14. 9p . (12.14[) . and the 
fact that u = for a = 0, to obtain that ■u(s) = iX[ 2 (s), which by (14. lip leads to 

"« - S ( ai'(,)»'-L(.)« ) = (*■") 

Its graph is shown in Figure EJ 

One can verify from the explicitly known asymptotic formulas for Ai that 

u(s) = —= - + 0(s~ 7/2 ) as s -> +oo. (4.15) 

On the negative real axis, u has an infinite number of zeros. These are the zeros of the 
Airy function Ai, and an infinite number of additional zeros that interlace with the zeros 
of Ai. 

Equations (14.91) and (I4.14p constitute the Schlesinger and (induced by it) Backlund 
transformations, respectively, for the case of Painleve XXXIV and applied to its zero 
vacuum solution (for the general theory of Schlesinger transformations see pTj; see also 
[21 Chapter 6]). 



4.3 Asymptotic characterization of the Painleve function u(s) 

We finally want to characterize the solution u of the Painleve XXXIV equation by its 
asymptotic properties. Recall that u is connected to the solution of the Painleve II 
equation q" = sq + 2q 3 — v with v = 2a + 1/2 by the formulas 

u(s) = 2- 1 / 3 £/(-2 1 / 3 s), U{s) = q 2 (s) + q'{s) + |. (4.16) 
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Figure 6: The solution of the Painleve XXXIV equation for a = 1. 

Assume that v > —1/2. It is shown in [31] (see also [2H Chapters 5, 11]) that the 
solution q(s) of the Painleve II equation corresponding to the Stokes multipliers (13.141) 
exhibits the following asymptotic behavior in the sector args £ > jF] 



/ ["+1/2] 

J-- b n {-s)-* n l 2 + 0(s-^ +l l 2 V 2 - 1 ) 

» 2 n=0 

+ c + (- S )-^^e-^ ? ^ 3/2 (l + 0(^ 1 / 4 )) 

2tt 4t 
~3~' ~3 

I [f+l/2] 

\K E &n(-«)- 3rl/2 + 0( S - 3[ ^ 1/2]/2 ) 

V 2 



g(sj 



or 



q(s) 



I 2tT 47T-, . . , 7T 7T-, 

as .s — oc . arg s £ ( — , —J , arg(-s) £ (--, -J . 



(4.17) 



n=0 



-s) 2 4 e 



^(-^(i + o^- 1 / 4 )), 



-27T 47T 



as s -> oo, args£[^,^), arg(-s) £ [--, -) 



3 ' 3 



3' 3' 



(4.18) 



where we have used the notation [r] for the integer part of the positive number r, i.e. 

[r] £ N , [r] < r < [r] + 1. 

The coefficients c+ and c_ of the exponential terms, which oscillate on the respective 
boundaries of the sector are given by the formulae 



7T 



and 



-7T(>+i)i 



7T 



■ +1 -rHr8 + „) 



(4.19) 



(4.20) 
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where T denotes the Gamma function. Moreover, either the relations (j4.17p . (14.191) or the 
relations (14.181) . (I4.20p can be taken as a characterization of the solution q(s). 

Alternatively, the solution q(s) can be characterized by its comparison to the Boutroux 
tri-tronquee solution q( tn - tron <i) ^ Q f the Painlevee II equation, which is defined as the 
unique solution satisfying the asymptotic condition 




CO 



s 

n=0 



4:71 

as s — > oo, arg s = n + arg(— s) G (0, — ) . (4.2T 

3 



The solution q(s) we are working with is the one whose asymptotic behavior as s — > — oo 
is given by the equation 



q(s) - q ( tri - tron i\ s ) = — — 2-?"-ira + u) 

71 

x | s |-I-? e -¥l«l 3/2 (1 + 0{s- 1 ^)) , as s - -oo. (4.22) 

The coefficients b n of the asymptotic series in (I4.17p . f!4. 18j) . and (14.211) are determined 
by substitution into the Painleve II equation. Indeed, the following recurrence relation 
takes place 

Qn2 _ 1 n+1 1 n+1 n+2-l (4.23) 

X- 

m=l 1=1 m=l 



n+1 ^ n+1 n+2-L 

';;+2 ~ b n — 2^ b m b n+2 -m ~ ~ °ib m b n+ 2-l-m- 



Using relation (14.161) between the Painleve II and Painleve XXXIV functions we arrive 
at the asymptotic characterization of the function u(s) of Theorem II. 51 

Proposition 4.1 The solution u(s) of the Painleve XXXIV equation which appears in 
Theorem \1.5[ is uniquely characterized by one of the following asymptotic conditions 



[2a+l] 

+ d +s ~ 3a+ h-l s3/2 {i + o(s- 1 / 4 )), 



MS 



, 7T 7T-i 

as s -> oo, args G (--, -J, (4.24) 



or 



u{s) 



[2a+l] 

° + ST anS - 3 -^ +0 ( a -3[2a + i]/2-n 
+ C /_ S - 3a+ 5e-! s3/2 (l + 0( S - 1 / 4 )), 



as s -> oo, args G [-^, ^), (4.25) 
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where 



d± = -2~ 6a -^T(l + 2a). (4.26) 



7T 



Alternatively, the solution u(s) can be characterized by the asymptotic relation 



-2am 



U[S 



71 



X s 



- 3Q +i e -| s 3/^ 1 + ^-i/4)^ ass ^ +00 . (4.27) 



The Painleve XXXIV tri-tronquee solution determined by the asymptotic 

condition 



u (tri-tronq) / \ ^ 



oo 

Q! ^ — V 3n + l _ / 71 "\ /, no i 

+ > a n s 2 as s —> oo, arg s G ( — vr, — J . (4.28) 

n=l 



Finally, the coefficients a n of the asymptotic series above can be expressed in terms of the 
coefficients b n defined in ( 14. 23[) . with v replaced by 2a + 1/2: 

n+l 372 — 2 1 \ 

2 2 a n = o n+1 ^7^~ & ™ + 2 bkbm - 

k,m>l;k+m=n+l 

Remark 4.2 The leading asymptotics of the Painleve II function q(s) as s ^ +oo is 
known (see [32]; see also [2U Chapter 10]). Unfortunately, the leading term is not enough 
to derive the corresponding asymptotics as s — > — oo of the Painleve XXXIV function 
u(s). Indeed, the leading asymptotics of q(s) as s — >• +oo is of the form 

q(s)r,^ co th^s^ + X (4.29) 

(the phase x is known) and it cancels out in the right-hand side of equation (13.221) . The 
better way to study the large negative s asymptotics of the function u(s) is via the direct 
analysis of the model RH problem for \I/ Q . The case a = 1 shows that we might expect 
oscillating behavior as s — > — oo (see Figure E]) and indeed, assuming that a — 1/2 ^ No, 
we are able to show that 

u(s) = -^=cos \\{-sf /2 -an j + 0(l/s 2 ), as s -> -oo. (4.30) 



The proof of (14.301) will be given in a future publication. Moreover, we conjecture that 
asymptotics (14.301) determines the solution u(s) uniquely. 
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